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SYMBOLS  AND  NOTATION 

a  s  speed  of  sound 
A  s  through-flow  area 
B  =  dimensionless  parameter 
FX  =  con^nressor  forces  acting  on  the  fluid 
H  =  total  enthalpy 

L  =  effective  length  of  equivalent  duct 
m  s  mass  flow  rate 
NBM  =  dimensionless  parameter 
P  -  static  pressure 
PR  =  pressure  ratio  (P^  /  P,,) 

Q  -  rate  of  heat  addition 
t  -  time 

T  -  static  temperature 
TR  -  temperature  ratio  (T^  /  T„) 

U  -  mean  rotor  velocity 
u  =  through-flow  velocity 
V  =  volume 

W,*  =  rate  of  shaft  work 
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SYMBOLS  AND  NOTATION  (continued) 


Greek  Symbols 

p  =  density 

Xp  s  cotqnressor  time  constant;  lag 


V 
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flow  coefficient 
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Actual  Corrected  Speed 


Subscripts 

b  =  stage  bleed 
c  -  compressor 
d  =  discharge  duct 
i  -  combustor  flller  volume 
p  s  plenum 
ss  -  steady-state 
t  -  total  or  stagnation 


1.0  INTRODUCTION 


>  Paramount  to  the  previous  and  continuing  development  of  the  aircraft  gas  turbine 

engine  is  an  uiklerstanding  of  compressor  instability.  Cintiptessor  instabilities,  as 
evidenced  in  the  axial  flow  compressors  used  in  modem  jet  engines,  have  a  detrimental 
effect  upon  the  operating  behavior  of  such  engines.  Reduced  performance  and  engine  life 
are  only  two  detriments  of  compressor  instability,  but  are  considered  among  the  most 
important  ones  to  aircraft  engine  designers.  The  lack  of  understanding  and  predictability 
of  unstable  compressor  behavior  in  the  compression  systems  of  jet  engines  has  led  to 
inconvenience  and  even  disaster.  As  a  result,  research  continues  in  an  effort  to  better 
understand  the  instability  behavior  of  high-speed  axial-flow  con^ressors. 

This  report  was  written  with  two  separate  goals  in  mind.  The  first  of  these  was  to 
present  some  of  the  surge  data  obtained  from  a  high-speed  10-stage  axial  flow 
compressor.  By  using  some  fundamental  data  analysis  techniques  as  well  as  conqruter 
based  flow  modeling,  the  basic  surge  behavior  of  this  corrqrressor  was  documented.  The 
second  and  more  challenging  goal  of  this  research  was  to  discover  the  reason  for  the 
strange  nature  of  the  surge/rotating  stall  boundary  in  this  conqnessor.  Changes  in  the 
discharge  plenum  volume  seem  to  have  had  no  effect  on  this  boundary,  an  observation 
which  runs  contrary  to  accepted  theory.  This  report  attempts  to  resolve  this  apparent 
contradiction  by  showing  that  certain  unexpected  flow  conditicMis  existed  during  surge 
which  may  have  affected  the  surge/rotating  stall  boundary  behavior  of  this  compressor. 
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1.1  Compressor  Behavior  and  Instability 


The  puipose  of  the  con^ressor  in  a  modem  jet  engine  is  to  increase  the  total  pressure 
entering  the  ctmnbostion  chamber  so  that  the  working  fluid  (air)  can  be  more  readily  and 
efficiently  heated.  Two  princ^al  types  of  compression  systems  are  typically  used  to 
perfonn  this  task.  Centrifugal  compressors,  which  utilize  centrifugal  force  to  help 
compress  the  working  fluid,  are  usually  characterized  by  a  relatively  high  pressure  rise 
per  individual  stage.  Axial  compressors,  which  rely  on  the  pressure  rise  produced  across 
rotating  airfoils  to  compress  the  working  fluid,  are  usually  characterized  by  high  mass 
flow  rates.  Generally,  centrifugal  corr^ressors ‘produce  higher  pressure  ratios  per  stage 
than  do  axial  compressors,  but  at  the  cost  of  added  weight  and  increased  total  pressure 
loss  due  to  flow  turning  between  individual  stages.  Due  to  the  weight  limitations  and 
high  mass  flow  requirements  imposed  upon  present  day  engines,  axial  compressors  have 
become  the  predominant  type  of  coirpression  system  employed  in  modem  jet  engines. 
Thus,  from  a  aircraft  propulsion  stand-point,  die  stability  behavior  of  axial  flow 
compressors  is  of  great  interest. 

The  performance  of  all  compressors,  including  the  axial  type,  is  often  described  by 
what  are  called  a  performance  characteristics.  Performance  characteristics  relate 
compressor  operating  variables  such  as  efficiency,  pressure,  temperature,  mass  flow,  etc, 
to  one  another.  As  an  example,  consider  the  "typical"  compressor  performance 
characteristic  for  stalled  and  unstalled  operation  shown  in  Figure  1.1.  This  performance 
characteristic  relates  the  total  pressure  ratio  produced  across  die  compressor  to  the  mass 
flow  through  the  compressor.  The  resulting  curves  define  the  stable  operating  points  for 
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Figure  1 .1  Typicar  compressor  characteristic 
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a  given  conqne^r  at  one  particular  conq>iessor  rotor  speed.  There  are  an  infinite 

I 

number  of  curves  for  an  infinite  number  of  compressor  speeds,  so  any  given  compressor 
will  have  a  set  of  conqnessor  characteristics  diat  vary  with  conqnessor  ^)eed. 

Generally,  there  are  two  operating  regimes  for  a  given  compressor.  In  Figure  1.1,  the 
curve  that  Includes  points  d  and  a  is  called  the  unstaUed  portion  of  the  compressor 
characteristic.  The  design  point  for  this  specific  compressor  (if  it  is  to  operate  at  this 
compressor  speed)  exists  somewhere  on  this  curve.  The  other  part  of  the  characteristic 
is  typically  called  the  stalled  portion  of  the  conpressor  characteristic.  When  the 
ctmpressor  c^nirates  at  some  point  on  this  curve  the  compressor  blades  are  stalled, 
causing  an  overall  reduction  in  compressor  performance. 

As  suggested  above,  it  is  desirable  for  the  compressor  to  operate  on  the  unstalled 
portion  of  die  compressor  characteristic.  However,  under  certain  conditions,  the 
conpressor  may  be  forced  to  operate  on  the  stalled  portion.  Consider  again  the 
compressor  pressure  characteristic  in  Figure  1.1.  Let  us  suppose  diat  die  compressor  is 
currently  operating  on  the  unstalled  portion  of  the  characteristic  at  an  operating  point 
somewhere  near  d.  Then,  holding  the  speed  constant,  let’s  further  suppose  that  the  mass 
flow  is  reduced.  In  a  conpressor  test  situation  where  the  compressor  is  not  part  of  an 
engine,  diis  is  usually  acconplished  with  a  valve  or  throttle  downstream  of  the 
compressor  exit.  In  this  case  the  tfarotde  is  slowly  closed. 

As  the  throtde  is  closed,  die  compressor  operating  point  begins  to  move  up  the 
unstalled  portion  of  the  compressor  characteristic  towards  a.  Let  us  suppose  throttle 
closure  is  stopped  once  the  compressor  operating  point  reaches  point  a.  The 
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ccunpressor  will  continue  to  operate  steadily  at  this  point.  Notice  the  dashed  line  shown 
on  the  "typical"  ccnnpiessor  characteristic.  This  line  is  called  a  throttle  line,  and  like  the 
conq»essor  characterise,  it  defines  all  the  »able  operating  points  for  the  throttle,  but  for 
a  given  throttle  flow  area  rather  than  for  a  givoi  conqnessor  qpeed.  As  is  the  case  for 
cominessor  pressure  characteristics,  there  are  an  infinite  number  of  throttle  lines  defined 
by  an  infinite  number  of  throttle  areas,  but  let  us  concern  ourselves  with  the  throttle  line 
that  passes  through  point  a  only. 

If  we  close  die  throttle  an  infinitesimal  amount,  the  compressor  operating  point  is 
forced  off  of  the  unstalled  portion  of  the  characteristic.  If  compressor  behavior  was  such 
that  the  operating  point  could  instandy  change,  the  compressor  operating  point  would 
"follow"  the  throtde  line  to  the  point  where  it  intersects  the  stalled  portion  of  the 
compressor  characteristic.  The  compressor  would  then  operate  quasi-steadily  at  this  point, 
shown  as  point  b  on  the  figure. 

Unfortunately,  compressor  operating  points  cannot  change  instantly.  Indeed,  there  is 
a  definite  time  lag  between  unstalled  and  stalled  operation.  Furthermore,  under  some 
conditions,  the  conqiressor  never  reaches  a  stable  operating  point  on  the  stalled 
characteristic.  Due  to  operating  conditions  that  are  not  yet  completely  understood  but  wiU 
be  briefly  discussed  later,  con:q)ressor  performance  beyond  the  unstalled  portion  of  the 
characteristic  is  often  defined  by  two  related  but  very  different  inst{d>ility  phenomena. 
These  ate  called  rotating  stall  and  surge.  Although  the  bulk  of  this  report  is  concerned 
with  surge  behavior,  these  two  types  of  unstable  compressor  behavior  are  too  interrelated 
to  discuss  one  exclusive  of  the  other,  so  both  are  discussed  herein. 
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1.2  Rotating  Stall  and  Surge 

Rntating  Stall 

Rotating  stall  manifests  itself  as  a  volume  or  cell  of  stalled  air  flow  diat  rotates  about 
die  connpessor  centeriine  at  some  fraction  of  die  ctmqvessor  rotor  qieed.  This  stall  cell 
rotates  in  the  same  direction  as  the  rotor  rotation  and  can  cause  a  significant  amount  of 
air  flow  blockage  in  the  compressor.  The  basic  mechanics  of  rotating  stall  inceptimi  are 
quite  simple  and  were  first  detailed  by  Emmons  et  al.  [1]  in  1955. 

Ccmsider  a  single  stage  compressor  rotor  modeled  as  a  linear  row  of  compressor 
blades  as  shown  in  Figure  1.2.  These  blades  are  at  a  high  angle  of  attack  and  are  moving 
in  the  direction  shown.  If  a  flow  disturbance  occurs  near  blade  ii  which  causes  an 
increased  angle  of  attack,  such  as  exhaust  gas  ingestion  or  off-design  inlet  flow,  the  flow 
on  the  suction  side  of  the  blade  separates.  Tim  causes  a  decrease  in  airflow  through  the 
blade  passage,  increasing  the  angle  of  attack  of  blade  i  and  decreasing  that  of  blade  iii. 
This  causes  blade  i  to  stall  and  at  the  same  time  inhibits  stalling  of  blade  iii.  After  blade 
1  stalls,  the  next  blade  in  the  blade  row  stalls  while  the  origiiuilly  stalled  blade  (blade  ii) 
tends  to  recover.  This  action  occurs  over  and  over  again  along  the  blade  row,  pn^agating 
in  die  direction  shown  in  the  figure.  If  the  initial  disturbance  is  strong  and  persists,  a 
fully  developed  rotating  stall  cell  may  develop. 

A  fully  developed  rotating  stall  cell  is  defined  as  an  area  of  extensive  blade  stalling 
that  rotates  around  the  conqiressor  aimulus.  Over  50  percent  of  the  armulus  area  can  be 
taken  iq>  by  a  rotating  stall  cell.  Moreover,  multiple  stall  cells  are  possible.  Air  flow 
through  a  rotating  stall  cell  is  nearly  zero,  and  is  sometimes  even  reversed. 
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accounting  for  the  blockage  mentioned  earlier.  Furfoermote,  die  blade  and  case 
teaqieratures  can  become  very  high,  so  prolonged  operation  in  rotating  stall  can  serioudy 
damage  a  compressor. 

When  an  instability  occurs  diat  initiates  rotating  stall,  the  compressor  can  go  into 
eidier  fnogressive  or  abrupt  rotating  stall.  Progressive  stall  occurs  when  the  (xmqnessor 
operates  on  a  continuous  compressor  characteristic.  That  is,  when  there  is  no  abrupt 
boundary  between  stalled  and  unstalled  performance.  Abrupt  stall  occurs  when  the 
conqnessor  characteristic  is  discontinuous  and  there  is  an  abrupt  boundary  between  stalled 
and  unstalled  performance. 

The  "typical"  con^ressor  characteristic  shown  in  Figure  1.1  is  conducive  to  abrapt 
stall,  so  in  this  case  the  compressor  does  not  typically  transition  smoothly  into  stalled 
operation.  Rather,  the  operating  point  of  the  compressor  teixls  to  oscillate  around  a  point 
defined  by  the  intersection  of  the  stalled  characteristic  widi  throtde  characteristic. 
Eventually,  the  operating  point  stabilizes  at  that  intersectitm.  Figure  1.3  shows  a  "typical" 
rotating  stall  inception,  where  the  compressor  operating  point  "^irals"  in  toward  a  quasi¬ 
steady  (operating  point  on  the  throttle/stall  characteristic.  Hie  term  quasi-steacfy  is  used 
to  define  rotating  stall  behavior  because  although  the  average  compressor  performance  is 
oftoi  defined  by  a  single  operating  point,  the  stall  cell  is  still  rotating  about  the 
compressor  rotor  causing  significant  fluctuations  in  pressure,  mass  flow,  and  ten^rature 
along  the  circumference.  Thus,  rotating  stall  is  stable  cmly  in  a  one-dimensirHial  sense. 
Surge 

The  otiier  possible  compressor  instability,  surge,  is  characterized  by  large  fnessure  and 
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MASS  FLOW 


Figure  Rotating  stall 


mass  flow  fluctuatitnis  which  tend  to  be  cyclic  in  nature.  During  surge,  conq>iessor 
operatitni  follows  a  cyclic  path  similar  to  that  shown  in  Figures  1.4  and  l.S.  Mass  flow 
can  drop  significantly  (Figure  1.4)  or  even  reverse  (Figure  l.S).  Because  of  its  low 
frequency  and  ifaythmic  behavior,  conq>ressor  surge  that  occurs  in  an  aircraft  engine  is 
often  referred  to  as  "banging."  Indeed,  surge  can  be  a  very  violent  occurrence,  causing 
stmctural  damage  to  the  engine  as  well  as  reducing  perfoimance. 

Unlike  rotating  stall,  which  is  quasi-steady,  surge  is  an  unsteady  phenomena.  Surge 
tends  to  occur  at  higher  speeds  than  rotating  stall  in  multistage  compressors.  Further,  die 
frequency  of  the  surge  cycles  tends  to  be  an  order  of  magnitude  lower  than  the  rotatitmal 
frequency  associated  with  rotating  stall. 

Recnverahilitv 

Another  important  differoice  between  mge  and  rotating  stall  is  the  ability  of  the 
compressor  to  recover  from  each.  Experiments  have  shown  that  in  many  compressors, 
iiKluding  the  ten  stage  rig  [2],  surge  is  iinich  easier  to  recover  from  than  rotating  stall. 
If  the  cause  of  the  instability  is  removed,  a  machine  that  is  experiencing  surge  will  almost 
immediately  shift  back  to  steady,  unstalled  operation.  This  is  not  the  case  for  rotating 
stall,  where  removing  die  cause  of  the  initial  instability  is  often  insufficioit  to  shift  the 
conqitessor  back  into  unstalled  operation. 

On  the  compressor  characteristic  shown  in  Figure  1.3,  note  that  the  operating  path 
followed  by  the  conqiressor  as  it  enters  rotating  stall  never  crosses  the  unstalled 
characteristic.  However  in  Figure  1.4,  as  the  crnnpressor  surges,  the  operating  path 
crosses  the  unstalled  characteristic  every  cycle.  It  is  because  of  this  that  the  compressor 
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Figure  14  Surge 
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Rgurel.5  Surge,  reverse  flow 


can  more  easily  recover  finnn  surge  than  from  rotating  stall.  When  die  cause  of  the  wge 
is  ronoved,  the  operation  of  the  conqiiessor  c(»itinues  to  follow  the  surge  (grating  path 
until  it  crosses  die  unstalled  characteristic.  Since  the  initiating  disturbaiK^  no  Itmgec 
exists,  the  ctxiqitessor  recovers  and  begins  operating  on  the  unstalled  characteristic  once 
it  crosses  the  characterise. 

The  rotating  stall  operating  path  iMver  crosses  the  unstalled  characteristic,  so  the 
compressor  must  recover  by  different  means.  Usually,  this  is  accomplished  by  somehow 
increasing  the  mass  flow  so  that  the  conqiressor  operating  point  is  pushed  past  c,  the  near- 
recovery  point.  Recovery  can  also  occur  by  reducing  rotor  speed  and  restarting  the 
compressor.  These  sort  of  drastic  corrective  measures  are  even  more  necessary  if  the 
phenomenon  known  as  hysteresis  is  present.  In  Figures  1.3  through  1.5,  the  stalled 
characteristic  is  seen  to  actually  pass  under  and  ovedap  the  unstalled  characteristic.  This 
means  .that  in  the  ovetlsp  region,  there  are  two  possible  operating  points  for  a  given  mass 
flow.  This  overlap  is  one  symptom  of  hysteresis.  As  previously  mentioned,  one  way  of 
recovering  from  rotating  stall  is  to  increase  the  mass  flow.  By  doing  so,  the  compressor 
is  forced  off  the  stalled  characteristic  and  onto  the  unstalled  characteristic.  If  the 
compressor  exhibits  hysteresis,  die  operating  point  will  remain  on  the  stalled  characteristic 
past  point  a  in  the  figures,  which  means  the  mass  flow  must  be  increased  even  more  than 
it  was  just  before  it  stalled.  Hysteresis  associated  with  tte  ten-stage  rig  is  discussed  in 
greater  detail  by  Copenhaver  [2,3].  Suffice  to  say,  surge  is  a  mote  desirable  condition 
than  rotating  stall  from  a  recoverability  point  of  view. 
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2.0  REVIEW  OF  LITERATURE 


This  section  briefly  describes  some  of  the  highlights  of  past  research  into  con^nessor 
surge.  The  first  two  subsections  discuss  work  ccmceming  the  surge^otating  stall  boundary 
ami  compressor  modeling  work  related  to  surge.  The  last  subsection  does  not  specifically 
concern  surge,  but  covers  research  that  was  done  using  data  from  the  same  ten-stage 
compressor  described  in  this  report.  Reviewing  previous  re^arch  on  this  compressor  is 
important  since  this  report  is  closely  related  to  most  of  that  work. 

2.1  The  Surge/Rotating  Stall  Boundary 

Surge  has  been  recognized  as  a  problem  in  turbo-compression  systems  since  the 
beginning  of  jet  engine  develppmoit.  Among  the  more  important  studies  done  on  this 
I^ienomenon  was  one  conducted  by  Emmons,  Pearson,  and  Gram  [1].  Enunons,  et  al. 
Investigated  the  relationship  between  compressor  surge  behavior  and  the  Helmholtz 
resonance  of  the  con^ression  system.  A  model  created  to  verify  this  theory  predicted 
surge  behavior  similar  to  experimental  results,  with  predicted  surge  frequencies  of  tlw 
same  order  of  magnitude  as  those  actually  observed.  Furthermore,  Emmons,  et  al. 
documented  the  speed  boimdary  between  surge  and  rotating  stall  (the  surge/rotating  stall 
boundary),  and  characterized  both  as  related  but  different  phenomena. 

The  surge/rotating  stall  speed  boundary  was  further  explored  by  Greitzer  [4,5]  using 
data  obtained  from  a  three-stage,  low-speed  compressor  rig.  Greitzer  developed  a 
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nonlinear,  luroped-paramem  model  to  predict  transient  c(mq>ressor  response  and  used 
data  from  the  three-stage  rig  to  verify  the  model  estimates.  Figure  2.1  is  a 
representation  of  the  model,  which  used  actuator  disks  in  constant  area  pq)es  to  simulate 
the  behavior  of  the  compressor  and  throttle.  Since  it  was  impossible  to  maAematically 
predict  conqnressor  pressure  rise  as  a  function  of  mass  flow,  the  unstalled,  stalled,  and 
reverse  flow  pressure  characteristics  of  the  compressor  had  to  be  provitted  from 
experimental  data.  Similarly,  the  time  lag  between  the  onset  of  compressor  instability  and 
the  establishment  of  fully  developed  rotating  stall  could  not  be  explicitly  calculated. 
Therefore,  the  time  lag  was  estimated  from  the  number  of  rotor  revolutions  that  occurred 
between  initial  compressor  instability  and  the  development  of  quasi-steady  rotating  stall. 
This  estimated  lag  was  then  used  with  a  simple  first  order  transioit  response  model  to 
simulate  the  time  lag  in  compressor  response. 

While  studying  the  surge/rotating  stall  boundary  with  the  model  and  compressor  data, 
Greitzer  found  that  the  location  of  the  boundary  depended  heavily  on  the  ratio  of  pressure 
forces  to  inertial  forces  in  the  conpiession  system.  This  ratio,  labeled  B,  comes  from 


the  equation: 


which  resulted  from  nondimensionalizing  the  governing  equations.  This  "B  parameter" 
is  also  related  to  the  Helmholtz  resonator  frequency  of  the  plenum  volume,  which  was 
deliberately  incorporated  into  the  model.  The  model  predicted  that  the  compression 
system  would  exhibit  surge  for  higher  values  of  B,  but  would  go  into  rotating 
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Figure  2.1  Greilzer's  compressor  model 
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stall  for  lower  values.  The  data  from  the  three-stage  rig  agreed  well  widi  the  motfel 
predictions. 

Small  and  Lewis  [6]  used  Greitzer’s  model  to  estimate  the  instability  behavior  of  their 
three-stage  high-q>eed  compressor  rig.  Using  the  B  parameter,  they  attenq>ted  to  shift  the 
surgefrotadng  stall  boundary  in  such  a  way  that  the  conq>ressor  would  stall  at  100%  speed 
with  nominal  IGV  settings.  They  did  this  by  altering  the  exit  plenum  volume  in  the 
compressor.  The  equation  defining  the  B  parameter  shows  a  significant  relationship 
between  the  exit  plenum  volume  and  the  surge/rotating  stall  boundary.  However,  they 
found  that  the  surgeAotadng  stall  boundary  shifted  only  slightly  (a  few  percent)  as  the 
exit  plenum  volume  was  varied. 

Hosny,  Steenken,  and  others  [7,8]  also  used  the  B  parameter  to  shift  the  surge/ 
rotating  stall  boundary  in  the  ten-stage,  high-speed  Energy  Efficient  Engine  (E^)  core 
compressor.  Earlier  builds  of  this  compressor  exhibited  rotating  stall  at  speeds  below 
60%,  but  the  researchers  wanted  to  observe  such  stall  behavior  at  higher  speeds.  Using 
the  equation  for  the  B  parameter  as  a  guide,  they  lowered  the  volume  of  the  exit  plenum 
in  an  attempt  to  raise  the  compressor  speed  defining  the  surg^totating  stall  boundary. 
As  expected,  the  boundary  did  shift  to  a  higher  q)eed,  but  it  did  not  precisely  follow  B 
parameter  equation.  A  plot  produced  in  one  the  reports  by  Hosny  and  Steenken  [8] 
showed  that  the  surge/rotating  stall  boundary  did  not  shift  as  much  as  the  equation 
estiiiuited  it  should  have. 

The  work  performed  by  Small  et  al.  and  Hosny  et  al.  demtmstrated  that  while  the  B 
parameter  estimated  the  surge/rotating  stall  boimdary  behavior  of  low-speed  compressors 
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(such  as  Greitzer’s  rig)  rather  well,  it  did  not  do  so  well  in  e^imating  such  behavior  in 
high-speed  compressors.  In  an  attempt  to  remedy  this,  Moore  and  Greitzer  [9,10] 
developed  a  nonlinear,  third-order,  lumped-parameter  conq>ressor  model  which  allowed 
for  the  simultaneous  occurrence  of  surge  and  rotiUing  stall.  They  showed  diat  the  lengdi- 
to-diameter  ratio  of  the  conqtressor,  in  conjunction  with  the  B  parameter,  was  necessary 
to  predict  high-speed  compressor  behavior.  This  is  important  in  multistage  designs  which 
tend  to  operate  at  high  speeds  and  have  increased  compressor  lengths. 
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2.2  Stage-by-Stage  Compressor  Model 


While  Mooie  and  Gieitzer  were  working  with  their  lunq)ed  parameter  model, 
Sugiyama  [11]  was  also  attempting  to  model  instabilities  in  multistage,  high-speed 
compression  systems.  Sugiyama’s  model  was  a  stage-by-stage  engine  model  which 
predicted  instability  behavior  in  an  entire  engine  using  individual  stage  characteristics 
instead  of  an  overall  con:q}ressor  characteristic.  Sugiyama’s  model  was  limited  to  short 
time  periods  (about  80  milliseconds)  and  had  some  instability  problems,  but  was  in 
general  an  excellent  early  attempt. 

Later,  Davis  [12]  developed  an  improved  one-dimensional,  stage-by-stage  dynamic 
conquressor  model.  Davis’  model  was  a  clear  in^rovement  over  Sugiyama’s  model  in 
that  it  allowed  for  longer  simulation  times  and  had  a  more  stable  numerical  solver. 
Furthermore,  the  Davis  model  could  simulate  recovery  from  surge  and  rotating  stall. 
Since  the  Davis  model  was  limited  to  modeling  the  compressor  rather  than  the  entire 
engine,  the  results  of  the  model  could  be  compared  to  avaUable  compressor  test  data. 

Davis  validated  the  model  by  comparing  the  model  estimates  to  three-stage  compressor 
data  provided  by  Greitzer  and  Moore  [5],[10]  and  to  ten-stage  data  provided  by  Hosny 
et  al.  [7,8].  Unfortunately,  detailed  characteristics  for  individual  stages  in  the  E’ 
confessor  were  not  availtUile.  Therefore,  Davis  synthesized  characteristics  for  a  nine- 
stage  compressor  using  a  stage-stacking  model  for  the  unstalled  characteristics  and 
estimations  based  on  low-speed  compressor  data  he  had  available  for  the  stalled 
characteristics.  The  model  predictions  con^ared  well  with  the  data  qualitatively,  but 
quantitative  con:q)arison  was  not  possille. 
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2.3  Ten-Stage  Compressor  Research 

In  addition  to  diis  repoit,  a  great  deal  of  previous  research  has  been  done  using  data 
from  diis  ten-stage  conqxressor  test.  The  first  major  research  effort  based  on  this 
conqnessor  data  was  an  Air  Force  technical  report  by  Copenhaver  [2]  publidied  in  1990. 
Besides  describing  the  general  behavior  of  the  conqnessor  during  instability,  Ct^ienhaver 
developed  the  stall  and  unstalled  conq)ressor  stage  characteristics  that  are  still  being  used 
in  other  research  efforts.  Copenhaver  also  validated  the  stage-by-stage  model  developed 
by  Davis  [12]  using  the  10-stage  data. 

In-stall  hysteresis,  which  is  the  mass  flow  required  to  move  from  the  initial  in-stall 
operating  point  to  the  near-recovery  point,  was  investigated  by  Copenhaver  and  Okiishi 
[3].  This  conq>re8sor  possessed  an  unusually  large  degree  of  in-stall  hysteresis,  which 
was  related  to  the  stalling  problems  associated  with  the  engine  from  which  this 
compressor  design  was  based.  Choking  in  the  rear  stages  of  the  compressor  was 
theorized  to  be  one  cause  of  this  large  degree  of  hysteresis. 

Later,  Boyer  and  O’Brien  [14]  used  the  conqnessor  characteristics  compiled  by 
Copenhaver  [2]  to  further  e;q)lore  die  hysteresis  and  recoverability  of  this  ctnnpressor.  To 
this  end,  they  used  the  Davis  model  [12]  to  simulate  conqiressor  stall  and  recovery. 
Using  simuli^ted  Laterstage  bleeds  and  increases  in  flow  padi  area  to  relieve  rear-stage 
choking,  they  determined  that  the  recoverability  of  the  compressor  seemed  to  depend  on 
the  location  of  the  bleed  flow.  Their  efforts  also  fiudier  validated  Davis’  stage-by-stage 
model. 

Gorrell  [IS]  built  on  Copenhaver’s  work  by  continuing  to  explore  the  rotating  stall 
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behavior  of  this  conq»essor.  By  reviewing  the  discharge  plornm  data,  Goirell  was  able 
to  deteiniiiie  dre  annular  extent  of  the  rotating  stall  cells.  Significant  transition  regims 
and  even  areas  of  reversed  flow  were  evident  in  those  stall  cells.  More  recently,  Bloch 
[16]  used  dre  ten-stage  data  to  verify  a  stage  characteristic  prediction  mednxL  This 
method  used  blade  angles  and  flow  geometry  to  jnedict  the  stage  duyracteristics  of  a 
given  cominessor.  The  results  of  this  method,  as  iqjplied  to  the  ten-stage  compressor, 
were  promising. 


21 


3.0  TEST  COMPRESSOR  DESCRIPTION 


The  data  presented  in  diis  pqwr  were  obtained  from  a  compressor  test  carried  out  at 
die  Compressor  Research  Facility  (CRF)  at  Wright  Patterson  AFB,  Ohio.  The  test 
ctnnpressor  design  was  based  on  a  10-stage  axial-flow  compressor  core  from  a  modem 
aircrafr  gas  turbine  engine.  This  compressor  had  a  (ksign  pressure  ratio  of  8.3,  a  design 
corrected  speed  of  10913  ipm,  and  a  design  corrected  mass  flow  of  24.69  kg/sec  (54.44 
Ib/sec). 

The  compression  system  in  the  actual  aircrafr  engine  consisted  of  both  a  3-stage  low- 
pf^'^ssure  fan  and  a  10-stage  hi^-piessure  compressor  along  with  the  necessary  ducting. 
The  test  compression  section  built  from  the  engine  did  not  include  the  fan,  but  it  did 
include  the  inlet  case  section,  the  high-pressure  compressor,  the  diSiiser,  and  pan  of  the 
combustor  section.  This  test  compression  section  is  shown  in  Figure  3.1. 

Inlet  and  discharge  sections  were  designed  and  built  so  that  the  test  compression 
section  could  be  incorporated  into  the  CRF  test  cell.  The  resulting  test  compressor  is 
shown  in  Figure  3.2.  While  the  bypass  flow  tiiat  occurs  in  the  actual  engine  was  not 
simulated,  the  flow  path  was  designed  to  simulate  air  flow  through  the  actual  engine  as 
much  as  possible.  At  the  same  time,  tlw  ability  to  vary  botii  the  flow  path  and  the 
conpressor  geometry  was  made  a  requirement.  To  this  end,  the  inlet  guide  vanes  and  the 
stators  in  die  first  two  stages  could  be  moved  independently  to  provide  varUd>le  blade 
stagger  angles  and  thus  variable  flow  conditions.  Moreover,  provisions  for  internal  and 
external  bleeds  from  die  compressor  core  were  also  included  in  the  design.  It  is  noted 
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MTemtAI.  FLOW  LOSSES 


Figure  3.1  Test  compressor  compression  section 


Hgure  3.2  CRF 1 0-stage  test  compressor  diagram 
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here  dut  data  presented  in  this  report  were  obtained  with  all  vanes  and  stators  held  at  a 
constant,  ixmiinai  angle.  Also,  the  external  bleed  was  closed  off  and  the  internal  bleed 
was  open  in  all  of  die  data  presented  in  this  report. 

The  test  cotn|»essor  inlet  secdtm  is  shown  in  Figure  3.3.  The  inlet  section  was 
designed  so  duu  die  airflow  going  into  die  conqiression  section  would  be  similar  to  the 
airflow  going  into  the  high-pressure  compressicMi  secticm  of  an  actual  engiru;.  Preswirl 
vanes  helped  simulate  airflow  angles  that  occur  as  a  result  of  the  fan  in  the  actual  engirt. 

Actual  flow  path  simulation  and  variable  getnnetry  were  also  incorporated  into  the 
discharge  section.  Aside  fitmi  the  flll^  volume  and  the  discharge  valve,  which  simulated 
the  combustor  and  turbine,  respectively,  die  discharge  section  included  variable  volume 
hardware  whidi  allowed  for  testing  at  nominal  and  ±  20  percent  nominal  discharge 
volume.  A  view  of  the  discharge  sectitm  is  ^wn  in  Figure  3.4.  Unfortunately,  the 
variable  volume  hardware  made  it  inqxissible  to  simulate  actual  engine  airflow  through 
the  discharge  section,  but  diis  was  judged  a  fair  exchange  for  the  ability  to  change  the 
conqnessor  discharge  volume  during  the  test. 

Definition  of  the  individual  stages  in  this  test  con^ressor  was  done  in  a  manner  that 
goes  against  accepted  convention.  While  compressor  stages  ate  conventionally  defined 
as  a  rotor-stator  combination,  stages  in  this  test  compressor  were  defined  as  a  stator-rotor 
combinatioa  (see  Figure  3.5).  This  was  required  since  instrumentation  used  to  monitor 
interstage  performance  was  typically  mounted  on  the  leading  edge  of  the  stator  in  each 
stage. 

A  brief  description  of  the  buildings,  rooms,  computer  systems,  and  drive  systems 
ctmaprising  the  Compressor  Research  Facility  is  included  in  Apper^ix  A. 
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Figure  3.3  Test  compressor  inlet  section 
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Figure  34  Test  compressor  discharge  section 
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Figure  3.5  Test  compressor  stage  definition  and  blade  locations 
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STATOR  lO 


4.0  COMPRESSOR  INSTRUMENTATION 


In  order  to  document  ctmqnessor  performance  and  bdiavior,  the  test  con^ressor  was 
thoroughly  instrumented  to  i»ovide  pressures,  tenqwratures,  and  odier  performance  values 
during  testing.  Much  of  this  instrumentation  and  the  data  associated  with  it  was  not  used 
in  this  repon.  Therefore,  this  section  will  deal  only  with  the  instrumentation  used  to 
provide  data  for  this  report.  A  more  thorough  explanation  of  the  data  acquisition  methods 
used  during  the  test  is  presented  in  Appendix  B.  A  more  detailed  discussion  of  all 
instrumentation  is  provided  by  Copenhaver  [2]  and  Copenhaver  and  Worland  [17]. 

4.1  High-Response  Pressure  Measurements 

During  the  con:q}ressor  test,  the  data  acquisition  system  was  required  to  accurately 
record  pressure  fluctuations  of  50  Hz  or  less.  This  requirement  was  set  down  so  that 
rotating  stall  events,  with  fundamental  frequencies  of  around  SO  Hz,  could  be  documented. 
To  this  end,  high-response  pressure  measurements  were  acquired  with  frequency  responses 
of  200  Hz.  This  was  adequate  to  document  tlw  rotating  stall  events.  The  surge  events 
presented  in  this  report  occurred  at  frequencies  near  S  Hz,  so  the  200-Hz  frequency 
response  of  the  high-response  data  adequately  documented  the  surge  events  as  well. 

Transient  pressure  data  were  obtained  from  six  pressure  transducers  comprising  four 
individual  probes  located  in  the  compressor.  Three  of  these  probes  were  forward  and  aft 
facing  pressure  probes  called  Mach  probes.  The  other  reiraining  probe,  which  was 
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actually  a  siinple  pressure  tap,  measured  static  pressure.  These  four  probes  provided  all 
of  the  transioit  pressure  data  presented  in  this  report. 

Sketches  the  two  types  of  Mach  probes  used  in  the  test  are  shown  in  Figures  4.1  and 
4.2.  Each  Mach  probe  consisted  of  at  least  two  impact  tubes  and  one  thermocouple. 
High-response  pressure  measurements  were  obtained  through  the  in:q>act  tubes  by  pressure 
transducers  located  on  the  outside  of  the  ccnnpressor  case.  Forward-facing  and  rearward¬ 
facing  measurements  were  obtained  through  the  angled  impact  tubes.  The  angle  of  the 
impact  tubes  allowed  the  dynamic  (velocity)  pressure  component  of  the  airflow  to  be 
monitored.  Thus,  the  Mach  probes  were  configured  to  measure  the  total  pressine  of  the 
air  at  the  measurement  location.  The  value  of  having  forward-facing  and  rearward-facing 
total  pressure  measurements  will  be  demonstrated  later  in  this  report.  A  sketch  of  a 
static  pressure  "probe"  is  ^own  in  Figure  4.3.  As  was  the  case  with  the  Mach  probes, 
pressure  measurements  were  obtained  through  a  tube  by  a  high-response  transducer 
located  outside  die  conqiressor  case.  However,  the  tube  used  here  was  not  angled  as 
those  on  the  Mach  probes  were,  but  was  angled  perpendicular  to  tlie  air  flow.  In  this 
way,  only  the  static  pressure  was  measured. 

The  positioning  of  the  four  probes  in  the  test  compressor  is  shown  in  Figure  4.4. 
Further,  some  inqiortant  characteristics  of  each  probe  are  listed  in  Table  4.1.  Of  the 
seven  high-reqionse  pressure  transducer  outputs  available  firom  the  four  probes,  only  six 
were  used.  The  rear-facing  pressure  measurements  taken  from  the  inlet  Mach  probe  were 
found  to  be  erroneous,  so  data  provided  by  that  transducer  were  not  used.  All  high- 
response  data  were  stored  on  analog  tspe  during  the  compressor  test  for  future  recovery. 
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9*Staclc  Pressure 
B'Tocsl  Pressure 


High-response  pressure  probe  locations 


Table  4.1  High-response  Pressure  Transducers 


Pressure 

Transducer 

Location/Probe 

Transducer 

Type 

Pressure  Range 

Number  1 

Total  Press. 

Comp.  Inlet, 
Forward-Facing 
Mach  Probe 

Kulite  Model# 
XCQ-25 

172  kPa 

(25  psi) 

Number  2 
Static  Press. 

Stg.  10  Inlet, 
O.D.  WaU  Tap 

Kulite  Model# 
XTEL-190 

690  kPa 
(1(X)  psi) 

Number  3 

Total  Press. 

Comp.Discharge 
Forward-Facing 
Mach  Probe 

Kulite  Model# 
XCQ-150 

1030  kPa 
(150  psi) 

Number  4 
Total  Press. 

Comp.Discharge 
Rearwd-Facing 
Mach  Probe 

Kulite  Model# 
XCQ-150 

1030  kPa 
(150  psi) 

Number  5 
Total  Press. 

Comp.Discharge 
Forward-Facing 
Mach  Probe 

Kulite  Model# 
XCQ-150 

1030  kPa 
(150  psi) 

Number  6 
Total  Press. 

Comp.Discharge 
Rearwd-Facing 
Mach  Probe 

Kulite  Model# 
XCQ-150 

1030  kPa 
(150  psi) 

_ — _ 1 
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4.2  Time-averaged  Measurements 

Accurate  steady-state  measurements  were  required  during  the  cmiqiressor  test  in  order 
to  facilitate  several  goals  of  the  test  program.  One  of  these  goals  was  the  vdltage-to- 
engmeeting  unit  conversion  of  the  high-re^xmse  transducer  measurements,  discussed  in 
detail  in  Section  5  of  this  report.  Time-averaged  data  acquisitirm  was  used  to  obtain  the 
steady-state  measurements  needed  to  achieve  these  goals.  Hie  measuremoit  uncertainty 
of  the  time-averaged  data  was  much  inqiroved  over  the  high-re^>onse  data  uncertainty 
because  of  the  very  accurate  calibration  standards  used  and  the  large  number  of  sanqiles 
diat  were  averaged  to  make  each  measuremem  (see  Appendix  B).  Time-averaged  data 
acquisition  was  used  to  obtain  the  steady-state  measurements  because  of  the  lesser  (fegree 
of  data  uncertainty  associated  with  it. 

Sixty'^iine  pressure  probes  and  thirty  total  temperature  probes  were  used  to  collect  the 
time-averaged  data  used  in  this  report.  All  of  these  were  placed  at  either  the  conqitessor 
inlet,  the  stage  10  inlet,  or  near  the  two  discharge  locations.  Figures  4.5  and  4.6  show 
the  time-averaged  pressure  and  terrqierature  probe  locations  in  the  test  conqiressor.  The 
locations  of  the  higb-response  pressure  probes  are  also  shown  in  these  figures  for 
reference  and  can  be  related  to  the  probe  positions  shown  in  Figure  4.4. 

Total  pressure  and  temperature  measurements  at  the  compressor  inlet  and  discharge 
were  takra  by  probes  similar  to  the  Mach  probes  used  to  obtain  the  high-response 
pressure  measurements.  Drawings  of  the  inlet  and  discharge  probes  used  by  the  time- 
averaged  acquisition  system  ate  shown  in  Figure  4.7.  These  probes  allowed  simultaneous 
measurements  to  be  taken  at  five  specific  radial  positions.  Each  probe  was 
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Figure  4.5  Time-averaged  instrumentation  locations 
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Figure  4.6  Time-averaged  instrumentation  locations 


posttioiied  on  the  outer  dumeter  of  the  compressor  case.  Aldiou^  the  designs 
were  identical,  seiMuate  probes  were  used  to  obtain  the  total  pressure  and  die  total 
temperature  measurements. 

Time<averaged  total  jnessure  and  temperature  measurements  at  the  inlet  to  stage  10 
were  taken  widi  probes  mounted  on  the  leading  edge  of  the  ninth  stator  ring.  For  eadi 
set  of  radial  measurements,  two  adjacent  stators  were  used  to  mmmt  the  five  different 
{vobes.  Two  locations  on  one  stator  and  three  on  the  other  were  instrumented  for  each 
set  of  measurements  (see  Figure  4.5). 

The  probes  used  to  obtain  die  dme-averaged  static  pressure  measurements  were 
vutually  identical  at  all  three  measurement  stations  used  in  the  test  ctnnpressor.  All  time- 
averaged  static  pressure  measurements  were  taken  with  probes  mounted  <ni  the  inner  and 
outer  diameters  of  the  test  compressor  flow-path.  The  time-averaged  and  high-response 
static  pressure  probes  have  similar  configurations,  so  except  for  the  transducer  the 
schematic  drawing  shown  in  Figure  4.3  is  teprerentative  of  both  probe  types. 

One  of  the  time-averaged  total  pressure  measurements  was  obtained  with  the  inlet 
Mach  probe.  The  schematic  of  the  inlet  Mach  probe  shown  in  Figure  4.1  shows  the 
second  pair  of  impact  tubes  used  to  obtain  two  time-averaged  total  pressure 
measurements.  As  widi  the  inlet  high-re^nse  pressure  measurement,  only  the  forward- 
facing  total  pressure  measuremem  was  used  in  this  report. 

All  time-averaged  pressure  measurements  were  taken  with  Scanivalve  model  ZOC-14 
(^to-£perate-£alibrate)  pressure  transducers.  TcHal  tenqierature  measurements  were 
taken  with  bead  type  thermocoi^les  constracted  frcrni  0.152  mm  (0.006  in)  diameter 
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Figure  4.7  Time-averaged  total  pressure  and  temperature  probe  designs 


chiomel-aiuniel  thermocoi^le  wire.  Total  pressure  measuremoits  were  obtained  through 
the  inq>act  tubes  with  the  actual  transducers  mounted  on  die  outside  of  the  conqiressor 
case.  ten^rature  measurements  were  obtained  widi  die  thermocouples  placed  just 
inside  of  die  impact  tubes  tm  the  total  tenqieratute  probes.  All  time-averaged  data 
acquisition  was  limited  to  a  frequency  response  of  0.5  Hz.  Time-averaged  data  were 
stored  on  digital  tape  for  future  recovery. 

4.3  Other  Measurements 

Mass  flow  through  the  ctmqiressor  was  comrolled  by  the  discharge  valve  located  just 
downstream  of  the  variable  volume  hardware  (see  Section  3.0).  Monitoring  the  discharge 
valve  flow  area  was  acconqilished  by  the  placement  of  three  potentiometers  along  the 
perimeter  of  the  discharge  valve.  Changes  in  potentiometer  output  voltage  indicated 
changes  in  the  position  of  the  discharge  valve.  By  knowing  the  area  of  the  discharge 
valve  when  fully  open,  the  flow  area  could  be  estimated  based  on  output  voltage. 
Discharge  valve  position  data  were  recorded  on  analog  ttqie  with  the  high-response  data. 

The  compressor  flow  rate  was  measured  by  a  calibrated  48.3-cm  (19-inch)  throat- 
diameter  flow  venturi  placed  approximately  31  meters  (25  p^  diameters)  downstream 
of  the  conq)tessor  exit.  Only  time-averaged  data  were  available  from  the  venturi. 
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5.0  DATA  REDUCTION  AND  INSTRUMENT 
CAUBRATION 

Most  of  die  data  recorded  from  the  10-stage  cooqvessor  test  woe  not  relevant  to  this 
rqwit.  For  diis  leascm,  the  first  task  was  to  determine  which  documented  evoits  would 
be  of  interest.  When  these  were  determined,  the  data  recorded  during  these  events  were 
recovered.  The  high-response  data,  which  were  recorded  on  analog  tape  in  the  form  of 
measured  voltages,  were  of  primary  interest.  To  use  the  data,  each  recorded  data  set  had 
to  be  digitized  and  calibrated.  This  process  of  data  recovery,  reduction,  and  calibration 
is  described  in  the  following  three  subsections. 

5.1  Data  Recovery 

This  report  was  limited  to  the  presentation  and  discussion  of  compressor  surge 
phenomena.  Therefore,  only  surge  events  were  of  interest.  During  the  10-stage  test,  at 
least  24  instances  of  surge  were  observed  (see  Table  5.1).  Of  these,  data  were  available 
for  22.  Many  of  these  surges  occurred  when  the  flow  path  gecxnetry  was  off  nominal, 
such  as  different  stator  or  inlet  guide  vane  angles.  These  were  ignored  since  the  effects 
of  variable  flow  padi  geometry  on  surge  behavior  were  not  die  subject  of  this  report. 
After  removing  these  from  consideration,  there  remained  13  documented  and  recorded 
surge  events  at  a  variety  of  speeds  and  plenum  volumes.  One  other  surge  event  not 
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included  in  the  table  was  too  short  to  be  considered  useful  and  was  not  recovered'. 


Therefore,  tlw  final  result  was  12  useful  and  recoverable  surge  events.  AU  of  these 
occurred  at  either  81,  82  or  87  percent  corrected  speed. 


'Thb  evert  was  imerestlng  because  K  begin  as  a  surge  but  after  one  cyde  ended  up  In  rotating  staH.  In  the  past,  this  author 
has  noted  simiiar  behavior  In  a  rMferert  compressor,  where  the  discharge  valve  was  dosed  far  enough  to  drive  the  compressor 
past  surge  and  into  rotating  stal.  Rtepossiito  that  discharge  valve  dosme  during  the  lO-stage  test  caused  the  surge-rotating  stall 
transition  in  the  same  way. 
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5.2  Data  Reduction 


High-response  Hata  for  the  12  surge  events  were  recorded  by  an  FM  ti^  recorder  in 
analog  form.  In  order  to  use  this  data  on  cmiqwter  for  both  analytical  and  plotting 
purposes,  it  was  necessaiy  to  conveit  the  analog  data  into  digital  form.  This  was 
accon:q>lished  with  die  analog  ti^  digitizing  system  (ATDS)  built  at  the  CRF. 

The  heart  of  the  ATDS  was  the  Massconqi  555/02  con^uter  system.  The  computer 
controlled  or  communicated  with  all  of  the  subsidiary  hardware,  including  the  analog  t^ 
player  and  anti-aliasing  filters.  The  actual  analog-to-digital  conversion  was  performed  by 
the  computer  as  well.  Once  digitized,  all  data  were  written  to  conjurer  hard  disk  for 
storage. 

Conversion  of  data  from  analog  to  digital  form  requires  some  knowledge  of  the 
transient  nature  of  the  data.  If  there  are  pertinent  high  frequency  components  in  the  data, 
the  sample  rate  during  the  conversion  must  be  high  enough  to  resolve  them.  However, 
high  sample  rates  create  large  amounts  of  data  which  can  quickly  fill  all  available 
computer  memory.  Somewhere  between  these  two  requirements  is  an  optimal  sample 
rate. 

Generally,  it  is  desirable  to  have  at  least  four  sample  points  over  one  period  or  cycle. 
For  example,  if  the  fiequency  to  be  resolved  is  one  cycle  per  second  (1  Hertz),  the 
digitization  sample  rate  should  be  at  least  four  san:q>les  per  second.  It  was  decided  that 
200  Hz  would  be  the  maximum  frequency  of  interest  in  the  high-response  data,  so  a 
sample  rate  of  1000  sanq)le8  per  second  was  used.  This  provided  a  minimum  of  five 
samples  per  cycle  at  200  Hz. 
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OiKc  the  sample  rate  was  chosen,  an  anti-aliasing  cut-off  frequency  had  to  be 
detennined.  Aliasing  occurs  when  high  frequency  conqwnents  in  die  raw  analog  data  are 
digitized  and  appear  in  the  digitized  data  as  low  frequency  ctm^nents.  Aliasing  is  a 
pitfall  of  any  digitization  process,  and  is  avoided  only  by  using  a  sample  rate  at  least 
twice  as  great  as  the  highest  expected  frequency  component.  This  is  called  die  Nyquist 
rate.  By  filtering  out  any  frequency  components  higher  than  one-half  the  sample  rate, 
aliasing  is  avoided. 

Under  ideal  circumstances  an  anti-aliasing  cut-off  frequency  of  500  Hz  would  seem 
appropriate  in  this  case.  However,  the  nonlinear  behavior  of  the  low-pass  filters  used  to 
remove  the  higher  frequency  components  required  deeper  consideration.  The  low-pass 
filters  used  by  the  ATDS  attenuate  frequency  conqionents  higher  than  the  cut-off 
frequency  but  do  not  completely  remove  them.  If  a  cut-off  frequency  of  500  Hz  had  been 
used,  the  filters  would  have  allowed  frequency  components  slighdy  higher  than  500  Hz 
to  pass.  This  would  have  caused  aliasing.  Therefore,  a  cut-off  frequency  of  400  Hz  was 
used  since  this  guaranteed  a  significant  reduction  in  the  signal  strength  of  frequency 
components  higher  than  500  Hz. 

Finally,  1  of  the  12  analog  data  sets  was  digitized  again  at  a  higher  sanqile  rate.  This 
"super  digitized"  set  was  produced  to  provide  more  data  points  for  Fourier  transform 
analysis  of  the  rotating  stall  phenomena  observed  in  the  data  after  initial  review.  This 
data  set  was  digitized  using  a  sample  rate  of  10,000  samples  per  second.  All  other 
parameters  were  unchanged.  A  summary  of  the  digitization  parameters  is  given  in  Table 
5.2. 
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Table  S.1  Observed  and  Usable  Surge  Events 


VohHne 

Total 

hhanber 

of 

Surge 

Events 

Minus 

Surge 

Events 

wHhno 

data 

Minus 

Surge 

Events 

w/stator 

angle 

change 

Minus 

Surge 

Events 

w/IGV 

angle 

change 

Minus 

Surge 

Events 

w/dlsL 

screen 

remvd. 

Minus 

Surge 

Events 

w/bleed 

on 

Total 

Number 

of 

Usable 

Stage 

Events 

80% 

Nominal 

9 

9 

5 

3 

3 

3 

3 

Nominal 

14 

12 

12 

12 

10 

8 

8 

120% 

Nominal 

1 

1 

1 

1 

1 

1 

1 

Table  5.2  Digitization  Parameters 


Digitization 

Parameters 

Alltoelve 
Surge  Events 

One  Surge 
Event 

1  Sample 

I  Rato 

1000"^/^ 

10000*^/^ 

Maximum 
Frequency  of 
Interest 

200  Hz 

200  Hz 

Anti-ABasing 

Frequency 

400  Hz 

400  Hz 
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5.3  Calibration 


Before  the  high-ie^nse  instrumentation  was  installed  in  the  test  con^ressor,  bench 

*  _ 

calfl>rations  were  performed  on  all  the  high-response  transducers.  The  original  intention 
was  to  use  diis  calibration  data  to  convert  the  output  voltages  into  mgineering  units  after 
the  test.  However,  day-to-day  observations  demonstrated  that  the  transducer  output 
voltages  changed  under  the  same  known  pressures  over  several  weeks.  This  voltage  drift 
made  the  bench  calibrations  less  accurate.  After  reviewing  similar  data  from  this  test, 
Gorrell  [15]  found  that  drifts  of  approximately  7  kPa  ( ~1  psi )  for  the  total  pressures  and 
2  kPa  (  0.3  psi  )  for  the  static  pressures  had  occurred.  Because  of  this  voltage  drift 
phenomenon,  calibration  of  the  high-response  data  was  accomplished  by  using  time- 
averaged  data.  Bias  errors  associated  with  the  total  and  static  steady-state,  time-averaged 
measurements  were  ^rproximately  ±  3.5  kPa  (  -0.5  psi ),  which  were  smaller  than  the 
errors  associated  with  the  voltage  drift. 

During  steady-state  compressor  operation,  high-re^nse  and  time-averaged 
measurements  were  taken  periodically  and  simultaneously  to  establish  steady-state 
operating  points.  At  these  points,  both  time-averaged  engineering  unit  measurements  and 
high-response  voltage  measurements  were  recorded.  If  a  specific  time-averaged 
instrument  was  in  the  same  immediate  location  as  a  high-response  instrument  measuring 
the  same  property,  the  two  instrument  readings  could  be  related.  This  is  what  was  done 
to  calibrate  the  high-response  data. 

High-response  data  taken  in  the  compressor  discharge  could  not  be  identically 
matched  to  time-averaged  data  since  the  related  time-averaged  instrumentation  was  2.5 
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cm  ( ~1  in )  upstream  of  the  hi^-resp<»ise  instnimentaticm.  It  was  necessary  in  this  case 
to  estimate  the  conditions  at  the  high-re^nse  instrumentation  location  based  on  die  time- 
averaged  measurements  and  some  sort  of  flow  model.  Ismtn^ic  flow  was  assumed 
between  the  two  locations,  since  it  was  a  bladeless  duct  with  a  small  area  change,  and  die 
static  conditums  2.5  cm  downstream  of  the  time-averaged  instrumoitatiun  were  estimated 
assuming  isentropic  expansion.  Totel  pressure  was  unchanged  (isentrqiic)  between  the 
two  locations. 

Time-averaged  and  high-re^nse  data  taken  cm  the  same  day  as  the  occurrence  of 
specific  surge  events  were  found  and  chosen  as  calibration  points.  The  high-re^XHise 
voltage  data  were  digitized  at  1000  samples  per  second  for  500  milliseconds.  This 
produced  500  high-response  voltage  measurements  which  were  averaged  to  determine  a 
mean  voltage.  Maintaining  a  95%  confidence  interval  for  the  mean,  the  worst  case 
difierence  between  the  upper  and  lower  confidence  limits  was  3.5  kPa  (  -0.5  psi  ). 
Generally,  the  difference  was  on  the  order  of  0.7  kPa  (  -0.1  psi ). 

Linear  fits  using  the  method  of  least  squares  were  calculated  using  the  averaged  high- 
response  measurements  and  the  time-averaged  data  to  determine  calibration  equations. 
Usually,  three  or  four  time-averaged/high-xesponse  data  point  pairs  were  obtained  for  the 
process,  and  the  resulting  equations  came  very  close  to  matching  these  points.  Each  set 
of  equations  included  calibration  data  for  seven  different  high-response  measurements. 

Several  of  the  12  surge  events  occurred  on  the  same  days,  which  allowed  one  set  of 
calibration  equations  to  be  used  for  more  than  one  surge  event.  Each  of  the  12  surge 
events  occurred  on  one  of  5  test  days.  Furthermore,  the  voltage  drift  of  die  high-response 
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instnimoits  used  did  not  aceed  0.5%  over  any  particular  day.  Therefore,  ptecisitm  was 
<Mily  minimally  affected  if  the  same  calibration  data  were  used  to  calibrate  two  or  more 
different  sets  of  surge  data  vdiich  occurred  on  the  same  day.  A  total  of  5  sets  of 
calibration  equations  were  used  to  calibrate  all  12  surge  events. 

The  calibration  equatioiis  were  also  effected  by  precision  errors  associated  widi  the 
time-averaged  and  high-response  data.  An  exanq}le  of  this  effect  is  shown  in  Figure  5.1. 
In  this  figure,  calibration  lines  using  the  worst  case  precision  errors  are  plotted  to  foim 
a  calibration  error  band.  The  nominal  calibration  line  was  used  to  calibrate  voltage  data 
from  probe  number  three  on  die  sectmd  of  dw  five  test  days.  This  is  one  of  die  worst 
cases  out  of  all  the  calibration  data  used  in  this  report.  In  this  particular  case,  errors  on 
the  order  of  4  kPa  (  ~0.6  psi )  were  possible,  fri  general,  most  measurement  errors  were 
smaller  than  this.  For  a  mote  complete  discussion  of  the  precision  and  uncertainty 
calculations  associated  with  the  measurements,  the  reader  is  referred  to  Appendix  B. 
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Figure  5.1  Example  caibration  error  band 


6.0  Data  Analysis 


The  calibrated  surge  data  are  presented  in  this  section.  TaUe  6.1  lists  the  12  surge 
events  and  smne  of  Ifaeir  basic  diaracteristics.  The  7  sets  of  time  resolved  data  were 
plotted  for  each  of  the  12  surge  events.  A  cm^lete  coinpflatioo  of  these  plots  is  given 
in  Appoidix  C.  One  of  the  most  important  observations  made  concerning  this  data  is  the 
spparent  lack  of  a  relationship  between  the  surgeAotadng  stall  boundary  and  the  plenum 
volume.  The  possible  reasons  for  this  lack  of  a  relationdiq)  are  as  yet  unclear.  The 
surge/rotating  stall  boundary  is  foe  subject  of  an  investigation  detailed  in  Section  6.2. 

6.1  Overview 

A  single  cycle  firom  1  of  foe  12  documented  surges  is  foown  in  Figure  6.1.  It  was 
qrparent  from  this  figure  that  the  general  nature  of  surge  in  this  compressor  was  typical 
of  surge  behavior  observed  in  otho:  compressors.  An  instability  brought  on  by  some 
irregularity  m  normal  operation  of  the  ccm^nessor  resulted  in  foe  break-down  of 
conqnessor  ability  to  maintain  a  pressure  rise.  This  break-down  caused  a  rapid 
depressurization  of  foe  discharge  plenum  immediately  upstream  of  foe  compressor.  The 
compressor  pressure  ratio  fell  until  the  conqnessor  was  able  to  recover.  At  the  point  of 
recovery,  the  compressor  began  to  rq)ressurize  the  plenum  until  foe  pressure  reached  an 
amount  neatly  equal  to  the  original  value  before  the  plenum  depressurized.  Unless  the 
cause  of  the  instability  was  removed,  foe  plenum  depressurized  again.  This  cycle  of 
depressurization,  often  called  blowdown,  followed  by  repiessurization  repeated  itself  until 
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surge  cycle 


the  cause  of  the  instability  was  removed  In  Figure  6.1,  the  rapid  blowdown  of  die 
discharge  plenum  followed  by  a  slower  repressurization  period  is  cleariy  evident. 

Whoi  the  forward-facing  and  rearward-facing  pressure  probe  data  fircnn  the  disdiarge 
were  OMnpared  it  became  tqiparent  that  flow  reversal  had  occurred  for  brief  periods 
during  surge.  Figure  6.2  shows  forward  and  rearward  pressure  data  finmi  a  discharge 
Mach  probe  during  surge.  The  dynamic  pressure  component  was  perceived  only  by  the 
forward-facing  probe  during  steady-state  operation,  so  the  forward-facing  probe  normally 
measured  total  pressure  and  the  rearwatd-fudng  probe  normally  represented  static 
pressure.  This  is  shown  in  Hgure  6.2,  where  the  higher  total  pressure  (forward-facing) 
is  plotted  with  the  lower  static  pressure  (rearward-facing)  just  before  the  first  blowdown. 
However,  for  brief  instances  during  surge,  die  rearward-facing  probe  measured  a  higher 
pressure  than  the  forward-facing  one,  as  the  figure  shows.  During  these  periods,  the  total 
pressure  was  measured  by  the  rearward-facing  probe  and  the  static  pressure  by  the 
forward-facing  one  because  die  dynamic  pressure  con^nent  was  briefly  being  perceived 
by  the  rearward-facing  probe.  This  indicated  flow  reversal. 

Given  the  static  and  total  pressure,  it  is  possible  to  estimate  the  flow  Mach  number 
at  die  discharge  probe  location.  In  Figure  63,  the  flow  Mach  number  based  on  the 
average  forward  and  rearward  pressure  measurements  in  the  discharge  during  surge  (81% 
speed,  100%  volume)  is  shown.  The  negative  Mach  numbers  indicate  reverse  flow. 
Although  not  shown,  the  pressure  data  from  the  other  four  surge  events  at  the  same  speed 
and  volume  provided  for  Mach  numbers  of  similar  magnitude.  The  other  seven  sets  of 
surge  event  data  also  showed  evidence  of  reverse  flow,  so  the  fact  that  flow  reversal 
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igure  6.2  Forward-facing  and  aft-fadng  probe  pressure  measurements 
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occuned  is  supported  by  data  fimn  ail  12  surge  events. 

High-iesponse  mass  flow  data  were  not  available  for  this  research  effort,  so  plots  of 
pressure  as  a  fiinctitm  of  mass  flow  during  surge,  like  the  {dots  shown  in  Figures  1.4  and 
1.5,  could  not  be  created.  However,  since  the  discharge  flow  Mach  number  was  available 
from  die  pressure  data,  pressure  versus  Mach  number  plots  could  be  goierated.  One  such 
plot,  based  on  data  from  Figures  6.2  and  6.3,  is  shown  in  Figure  6.4.  By  relating  the 
Mach  number  to  the  mass  flow  and  comparing  this  figure  to  Figure  1.5,  the  movement 
of  the  compressor  operating  point  from  unstalled  to  stalled  operation  and  back  again  is 
cleady  evidmt  Althouj^  interesting,  this  type  of  plot  was  not  pardculariy  useful  in  this 
research  effort,  so  similar  plots  of  the  other  11  surge  events  were  not  generated. 

During  the  reverse  flow  periods,  what  appeared  to  be  rotating  stall  cells  were 
observed  in  all  12  surge  events.  One  of  these  is  identified  in  Figure  6.1.  The  annular 
length  of  these  cells  stretched  at  least  from  the  discharge  pressure  probes  to  the  inlet  to 
stage  10.  This  fact  is  demonstrated  in  Figure  6.5,  where  pressure  data  from  the 
compressor  discharge  are  compared  to  pressure  data  from  the  stage  10  inlet.  The  two 
data  collection  locations  are  nearly  in  line  with  one  another,  and  it  is  clear  that  the  cell 
passed  both  locations  at  almost  the  same  time.  The  other  static  pressure  probe  located 
in  the  discharge  also  recorded  die  cell  passage,  but  it  is  out  of  phase  with  the  other  two 
probes.  Using  die  surge  data,  an  attenqit  was  made  to  isolate  the  rotating  stall  pattern 
from  die  surge.  This  was  done  by  subtracting  a  representative  blowdown  from  actual 
blowdown  data.  The  representative  blowdown  was  created  by  obtaining  points  along  a 
blowdown  trace  recorded  by  the  pressure  probes  in  the  discharge.  The  points  were 
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Figure  6.5  Rotating  stall  phase  relationship:  stage  10  inlet  and  discharge  pressure  probes 


chosen  from  blowdown  regions  where  the  rotating  stall  was  absent,  so  that  the 
representative  blowdown  would  be  free  of  rotating  stall  components.  Since  pure  surge 
is  generally  an  axisymmetric  phenomenon,  the  points  obtained  from  probes  facing  the 
same  direction  were  combined  and  a  curve  fit  was  calculated.  The  curve  fit  was  used  as 
a  representative  blowdown. 

To  adjust  the  magnitude  of  the  static  and  total  pressure  readings  once  the  blowdown 
was  subtracted,  average  in-stall  pressure  measurements  taken  at  79.5%  compressor  speed 
were  used  as  reference  pressures.  Qearly,  using  these  pressures  as  references  is  not  an 
altogether  correct  practice,  but  in-stall  data  above  79.5%  does  not  exist.  For  this  reason, 
the  data  resulting  from  this  process  should  be  considered  qualitative  and  not  quantitative. 
The  results  of  this  process  using  one  set  of  probes  are  shown  in  Figure  6.6.  The  general 
shape  of  the  stall  cell  resembles  the  rotating  stall  plots  produced  by  Copenhaver  [2]  using 
lower  speed  data  from  this  compressor.  This  was  interesting  because  it  seemed  to  support 
the  premise  that  these  pressure  abnormalities  were  indeed  rotating  stall  cells. 

In  several  cases  the  stall  cell  rotated  about  the  compressor  twice  during  a  blowdown. 
Therefore,  sufficient  data  existed  to  estimate  the  rotating  stall  frequency  had  the 
compressor  not  recovered  from  surge.  The  estimated  stall  frequencies  for  each  surge  are 
shown  in  Table  6.1.  It  is  clear  that  the  stall  frequency  did  not  change  as  the  volume  was 
changed,  but  the  frequency  did  increase  as  the  compressor  speed  was  increased. 

All  12  suiges  were  induced  by  reducing  the  mass  flow,  accomplished  by  slowly 
closing  the  discharge  valve.  Since  surge  could  cause  significant  structural  damage,  the 
compressor  was  not  allowed  to  surge  for  very  long.  Therefore,  once  the  compressor  was 
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isolation  attempt 


Table  6.1  Surge  event  characteristics 


seen  to  be  surging,  the  discharge  valve  was  usually  opened  quickly.  Figure  6.7  shows 
the  time  history  of  the  discharge  valve  flow  area  during  a  surge  event.  The  line  is  only 
slightly  negatively  sloped  most  of  the  period,  until  the  point  where  the  discharge  valve 
was  rapidly  opened. 

Also  shown  in  Figure  6.7  is  the  time  history  of  die  discharge  pressure.  In  diis  figure, 
it  is  clear  that  the  discharge  valve  was  still  being  closed  after  surging  had  started.  It  was 
observed  that  the  time  length,  or  period,  between  surge  cycles  continued  to  decrease  as 
the  discharge  valve  was  being  closed.  This  behavior  is  clear  after  careful  examination  of 
Figure  6.7.  The  time  period  between  surge  cycles  decreased  because  the  plenum 
repressurized  faster  when  the  outflow  was  less.  This  becomes  clear  if  one  imagines  trying 
to  inflate  a  balloon  with  a  large  tear  in  it  versus  inflating  one  with  a  tiny  tear  instead. 

Although  this  behavior  was  not  unexpected,  it  was  important  to  realize  that  surge  was 
not  occurring  under  constant  conditions  because  of  this  phenomenon.  Surge  cycles  from 
the  beginning  of  a  surge  event  could  not  be  compared  to  cycles  at  the  end  on  an 
equivalent  basis.  This  fact  made  it  difficult  to  attach  significance  to  an  average  surge 
frequency  based  on  summing  all  of  the  surge  cycles  during  a  surge  event.  However,  it 
was  possible  to  conqrare  surge  data  from  the  beginning  of  one  surge  event  to  data  from 
the  begirming  of  another  surge  event  because  the  discharge  valve  area  that  initiated  surge 
was  faMy  consistoit  for  all  three  volumes  at  a  given  speed.  This  allowed  for  the 
calculation  of  surge  frequencies  based  on  the  periods  of  die  first  surge  cycles. 

Table  6.1  includes  the  surge  frequencies  calculated  in  two  different  ways.  The 
average  surge  frequency  is  calculated  by  averaging  the  periods  of  all  the  surge  cycles 
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vatve  time  history 


ihiring  the  surge  cvast.  The  tqTpareat  surge  frequoicy  is  an  estimate  based  on  the  first 
cycle  (Hily.  In  both  methods,  the  period  is  the  time  betwem  the  low  pressure  points,  or 
the  troughs.  It  was  found  that  detennining  the  time  period  trmigh-to-trougfa  was  much 
easier  and  more  consistent  than  tiying  to  detennine  die  peak-to-peak  period. 

After  reviewing  the  surge  frequency  data,  a  lelationsh^  between  the  suige  frequoicy 
and  the  plenum  volume  seemed  to  emerge.  Most  theories,  including  the  conqiiessor 
model  discussed  in  the  next  section,  predict  an  increase  in  surge  frequency  with 
decreasing  plenum  volume.  Indeed,  at  81%  speed,  the  surge  frequency  seemed  to 
increase  as  the  plenum  volume  decreased.  However,  at  82%  ^leed  the  opposite  trend 
seemed  to  occur:  as  the  volume  decreased  so  did  the  frequency.  This  apparent 
contradiction  is  probably  due  to  the  lack  of  available  data.  For  example,  at  81%  speed 
and  120%  volume,  the  one  surge  event  recorocd  had  only  four  complete  cycles.  Four 
cycles  is  insufficient  for  a  good  estimate  of  the  average  surge  frequency.  Furtiier,  the 
tq^aient  surge  frequency  at  this  speed  and  volume  is  based  on  only  one  cycle,  since  only 
one  instance  of  surge  occurred  under  those  conditions.  So  even  though  the  surge 
fiequency  probably  increased  with  decreasing  volume,  a  lack  of  data  made  it  difficult  to 
prove  this  conclusively. 

Pressure  data  fixnn  the  compressor  inlet  was  also  available  for  study.  Figure  6.8 
shows  a  typical  time  history  of  the  inlet  and  discharge  total  pressures  during  surge  as 
recorded  by  forward-facing  Mach  probes.  Unfortunately,  the  rearward-facing  Mach  probe 
at  the  inlet  worked  only  intermittently  and  did  not  provide  useful  data.  The  passage  of 
pressure  waves  through  die  inlet  during  each  surge  cycle  is  evident  by  the  pressure  spikes 
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ire  6.8  Inlet  and  discharge  pressure  comparison 


in  the  figure.  There  is  a  definite  time  lag  between  these  pressiie  spikes  and  beginning 
of  blowdown  registered  by  the  discharge  Mach  jvobe.  It  was  possible  to  use  this  time 
lag  in  ccmjunction  with  steady-state  tempoature  and  presaire  data  to  estimate  the 
q)pioximate  locatitm  in  tibe  compressor  where  dte  initial  stall  that  caused  the  surge 
occurred.  In  all  12  cases  (81%  to  87%  speed),  the  initial  stall  seoned  to  originate  in  the 
last  3  or  4  stages  of  the  compressoT.  This  agreed  weU  widi  eailier  research  efforts 
ctmceming  this  cottqrressor  including  Copenhaver’s  [2]  work,  which  showed  that  ^alls 
near  78.5%  con^ressor  speed  originated  in  the  seventh  stage. 

Further  study  of  lower  ^>eed  data  from  this  ctnnpressor  showed  that  78.5%  speed  was 
the  highest  speed  at  which  rotating  stalls  were  recorded  in  this  c(»iq>ressor^.  Under 
nominal  through-flow  conditions,  diis  maximum  speed  for  stall  did  not  change  when  die 
plenum  volume  was  varied.  Further,  it  was  q>parent  from  the  surge  data  presented  herein 
that  the  minimum  q)eed  for  surge,  80.4%  (-81%)  corrected  speed,  also  did  not  change 
when  the  plenum  volume  was  varied.  Thus,  the  speed  boundary  between  surge  and 
rotating  stall,  shown  in  Figure  6.9,  did  not  change  when  the  plenum  volume  was  increased 
or  decreased. 


*  This  does  not  inciude  the  surge-to-rotating  stal  transition  mentioned  earlier  in  this  report  (see  foohiote  ’). 
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Figure  6.9  Surge/rotating  stall  boundary 


6.2  Surge/Rotating  Stall  Boundary  Investigation 


All  c(xiq>ie8sion  systems  that  exhibit  both  surge  and  rmating  stall  over  portions  of 
iheir  q)erating  range  have  what  is  called  a  surg^rotadng  stall  boondaiy  (hencefoith  to  be 
called  the  S/RS  boundary).  This  is  typically  defined  as  a  ^)eed  boundary  that  separates 
the  ctnnpressor  stall  regions  where  either  surge  or  rotating  stall  occur.  Whem  a  given 
c(»ipressor  operates  beyond  the  stall  line  with  constant  geometry  and  discharge  volume, 
surge  occurs  at  conpressor  speeds  above  the  S/RS  boundary  and  rotating  stall  at  speeds 
below  it. 

The  compressor  speed  defining  the  S/RS  boundary  is  unique  to  a  paiticular 
conpressor.  Some  of  die  factors  that  deteimine  this  boundary  are  blade  geometry,  the 
size  and  length  of  the  conpressor,  and  the  number  of  stages,  but  the  factor  that  receives 
the  most  attention  herein  is  the  conpressor  discharge  volume.  Most  dynamic  compressor 
models,  including  the  ones  to  be  discussed  later  in  this  section,  show  a  relationship 
between  the  discharge  volume  and  the  S/RS  boundary.  This  has  much  to  do  with 
conservation  of  energy  within  the  system.  When  pressurized,  die  discharge  volume 
behaves  much  like  a  capacitor  in  an  electrical  system,  storing  mechanical  potential  energy 
in  the  form  of  pressurized  air.  If  the  volume  is  large,  the  amount  of  potential  energy 
stored  is  greater  dian  if  the  volume  were  smaller.  It  has  been  found  that  the  more  energy 
there  is  in  the  system,  the  greater  the  likelihood  that  surge  will  occur.  Therefore,  larger 
discharge  volumes  toid  to  cause  sui^ge  instead  of  rotating  stall. 

The  10-stage  conpressor  tested  at  the  CRF  exhibited  both  surge  and  rotating  stall 
during  unstable  operation.  As  expected,  rotating  stall  tended  to  occur  at  lower  conpressor 
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speeds  while  surge  occurred  at  higher  speeds.  However,  discharge  voIuiim  changes  did 
not  seem  to  affect  the  S/RS  boundary.  In  Section  6.1,  it  was  shown  that  the  boundary 
remained  near  80%  conected  rotor  speed  despite  discharge  volume  changes  of  ±  20% 
iKHniiud.  The  surge  frequency  should  also  have  changed,  but  as  mentioned  in  the 
previous  sectitm,  the  surge  frequency  changed  very  little  (if  at  all).  Possible  reasons  for 
this  lack  of  a  volume  -  S/RS  boundary  relationship  are  investigated  in  this  sectitm. 

In  the  process  of  investigating  the  S/RS  boundary  behavior,  a  dynamic  compressor 
model  is  used.  After  comparing  the  estimated  (or  simulated)  and  acmal  behavior,  a  lack 
of  agreement  between  the  two  is  tqppareiit.  The  possibility  that  more  volume  was 
apparent  to  the  compressor  than  previously  thought  is  explored  in  detail,  and  several 
theories,  which  attempt  to  account  for  this,  are  formulated.  Several  other  hypotheses 
which  might  explain  the  S/RS  boundary  behavior  are  presented  at  the  end  of  this  section, 
but  none  of  these  are  investigated  in  detail. 
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6.2.1  Pre-t98t  S/RS  boundary  prediction:  the  NBM  parameter 


Before  the  test,  die  builder  of  the  test  ccmipressor  (Pratt  &  Whitney)  made  an  estimate 
of  the  surge/totating  stall  boundary.  This  was  (kme  to  aid  in  the  design  of  the  variable 
volume  hardware.  Based  on  this  surgeAt^adng  stall  boundary  estimate,  the  maximum  test 
plenum  volume  was  determined.  The  maximum  ten^ierature  allowed  [18]  {xit  limitations 
on  the  minimum  volume,  but  the  boimdaty  estimate  helped  the  designers  choose  a  factor 
of  safety.  Although  70%  nominal  volume  was  possible  based  on  temperature  limitiuions, 
80%  nominal  volume  was  diought  to  be  con^atible  with  the  overall  test  objective. 

To  make  the  boundary  estimate,  the  manufacturer  used  a  dimensionless  parameter 
related  to  Greitzer’s  B  parameter  [4,5].  The  NBM  parameter  [19],  was  defined  in  the 
following  manner 


NBM* 


PRa  y, 
2«^  > 


where: 

PR  =  pressure  ratio  at  stall 
Ug  =  axial  flow  velocity, 
a  =  speed  of  sound 
Ag  =  conq)ressor  flow-through  area 
=  plenum  volume 

Lg  =  effective  length  of  equivalent  duct  for  compressor 


Using  this  relationship  in  conjunction  with  their  previous  experience  with  similar 
compressors,  P&W  made  the  surge/rotating  stall  boundary  estimate  drown  in  Figure  6.10. 
At  any  point  on  the  S/RS  boundary  prediction,  NBM  7.  Obviously,  the  prediction 
failed  to  agree  with  the  actual  boundary.  This  disagreement  was  the  first  indication  that 
the  S/RS  boundary  behavior  of  this  compressor  required  further  investigation. 
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Fiqure  6.10  Pre-test  S/RS  boundary  estimate 


6.2.2  Stage-by-stage  model 

During  the  eariy  1990’s,  a  stage-by-stage  post-stall  c(mq>iessor  model  was  developed 
at  Arnold  Engineering  and  Development  Center  (AEDC)‘based  m  the  work  of  Dr.  Milt 
Davis  of  Sverdnq)  Technology  Inc.  [20].  This  model  was  designed  to  deal  with  the 
multiple  stages  and  high  pressure  ratios  associated  with  noodem  cmnpiessors.  The 
Dynamic  Turbine  Engine  Coti^ressor  Code,  or  dyntecx:,  was  and  still  is  the  state-of-the- 
art  in  con^ter  based  dynamic  con^iessor  modeling. 

DYNTECC  models  the  compressor  as  a  series  of  control  volumes,  much  like  tlw 
configuration  shown  in  Figure  6.11.  Individual  con^ressor  stages  as  well  as  ducting 
segments  are  modeled  as  elemental  control  volumes.  Mass  flows,  pressures,  bleeds,  blade 
forces,  etc.,  ate  related  to  one  another  in  each  control  volume  through  conservation  of 
mass,  momentum,  and  energy  principles.  Pressure  and  temperature  characteristics  for 
each  stage  are  requited  for  the  model  to  determine  the  blade  forces  and  draft  work. 

To  account  for  the  transient  response  of  the  compressor,  die  blade  forces  are  lagged 
using  the  first-order  lag  equation; 

D  dt  » 

and  a  specified  time  constant,  Tq.  The  value  of  Xp  depended  upmi  what  portion  of  the 
compressor  characteristic  the  compressor  was  operating  on.  At  stall  incepticm  (post-stall) 
up  to  some  minimum  mass  flow  coefficimt,  the  model  would  use  one  value  for 
During  pre-recovery  (or  near-recovery  in  the  case  of  rotating  stall),  Xp  would  take  on  a 
difierent  value  from  some  minimum  to  some  maximum  mass  flow  coefficient  value. 
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Control  Vohme 


Figure  6.11  Control  volume  for  stage-by-stage  model  (dyntecc) 
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Multipliers  included  in  the  input  to  the  model  controlled  the  mass  flow  coefficient  ranges. 
For  the  post-stall  value,  the  multipliers  were  based  on  die  near-stall  point  on  dw 
unstalled  characteristic.  In  the  case  of  the  pre-recovery  Xp  value,  the  mult^liers  were 
based  on  the  near-recovery  point  on  the  stalled  characteristic.  The  value  of  the  maximum 
flow  coefficient  multiplier  for  the  pre-recovery  Xp  value  dictated  the  de^ee  of  hysteresis 
in  that  particular  stage.  In  the  case  of  these  model  simulations,  pre-recovery  lag  was 
applied  up  to  a  flow  coefficient  value  1.2  times  the  near-recovery  point  value  on  the 
stalled  characteristics  of  all  10  stages.  This  flow  coefficient  value  was  based  on  the 
observed  hysteresis  of  the  con^ressor. 

The  reverse  flow  charaaeristics  could  also  be  manipulated  in  this  model.  Since  no 
reverse  flow  characteristics  exist  for  this  compressor,  deciding  what  these  characteristics 
should  look  like  is  up  to  the  researcher.  For  this  report,  the  reverse  flow  characteristics 
estimated  by  Bloch  [16]  were  used.  The  multipliers  for  slope  and  offset  in  the  model 
input  reflect  these  estimates.  All  of  the  specific  values  for  Xp  and  the  various  multipliers 
used  in  this  DYNTECC  model  are  shown  in  Table  6.2. 

As  previously  mentioned,  this  model  uses  the  individual  stage  characteristics  to 
simulate  overall  compressor  behavior.  For  conq>ressor  speeds  where  specific 
characteristics  do  not  exist,  the  model  uses  an  interpolation  routine  to  estimate  the 
appropriate  characteristics.  This  routine  uses  the  characteristics  from  a  higher  speed  and 
a  lower  speed  to  make  the  estimate.  In  the  case  of  this  con:q)ressor,  characteristics  for 
78.5%  speed  and  82%  speed  were  made  available.  The  characteristics  for  82%  speed  are 
shown  in  Figure  6.12.  Since  most  of  the  transient  data  recovered  were  for  speeds  near 
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Table  6.2  dyntecc  input  values 


Stage-by-stage  Model  Inputs 

Value 

Variable 

Names 

Percent  design  mass  flow 

37.70/41.65 

PWCDES 

Percent  design  rotor  speed 

78.5  /  81 

SPEEKIPC 

to.  Post-stall 

0.012  sec 

TLAOHPCT 

to.  Pre-recovery 

0.033  sec 

TLAOHFCB 

Post-stall  Tq  control;  Initial 

1.0 

PHIHPCm 

Post-stall  tp  control;  Final 

0.120  sii..  1-3 

0.235 

PHIHPCTL 

Pre-iecoveiy  Tq  control;  Initial 

1.2 

PHIHPCBR 

Pte-recovciy  Xq  control;  Final 

0.0 

PHIHPCBL 

Reverse  flow  characteristic  slope  for 
each  stage  (Pressure  only) 

4^A1.8^.5 

HPCROTAT 

Reverse  flow  characteristic  offset 
multiplier  (Pressure  coefficient) 

1.0  SlifH  1.3 
0.65 

PDROraPC 

Reverse  flow  characteristic  offset 
mult^lier  (Tenq)erature  coefficient) 

1.0 

TDROPHPC 

81%,  the  model  estimated  the  characteristics  at  that  speed  when  craipariscHis  betweoi  the 
model  and  actual  data  were  made.  actual  characteristics  used  were  con^iled  by 
Boyer  and  O'Brioi  [14]  and  Copenhaver  [2]. 

Most  of  the  work  involved  in  setting-iq>  the  10-stage  amqvessor  model  using 
DYNlECC  had  already  been  accmnplished  by  previous  researchers,  but  it  was  still 
necessary  to  adapt  the  model  to  this  research  effort.  Parts  of  the  con^ter  code  were 
altered  to  modify  the  output  for  tiM  plotting  software  used,  but  no  significant 
modifications  were  made  to  DYNTECC.  Once  set-up,  DYNTBCC  was  used  to  help 
accomplish  two  research  goals.  The  first  of  these  was  to  determine  if  the  model  could 
be  made  to  simulate  die  behavior  of  the  test  compressor  using  the  actual  con^ressor 
geometry  and  operating  conditions  as  irqiuts.  If  this  could  not  be  accomplished,  the 
second  goal  was  to  find  a  discharge  volume  which  would  allow  the  model  to  simulate 
actual  conqiressor  behavior  with  all  other  operating  conditions  held  constant.  The  first 
goal  was  not  met:  the  model  did  not  simulate  actual  behavior  using  nominal  input  data. 
However,  the  second  goal  was  accomplished.  The  efforts  made  and  results  from  both 
simulation  attempts  ate  described  next. 
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The  model  was  calibrated  at  100%  nominal  discharge  volunM  so  that  the  S/RS 
boundary  s^^ated  between  79%  and  81%  compressor  speed.  The  ctmpressor  time 
constants,  were  instrumental  in  acctnnplishing  this.  By  trial  and  error,  values  for  Xp 
were  found  which  caused  the  simulation  to  stall  at  79%  speed  and  surge  at  81%  speed 
with  nominal  discharge  volume.  These  were  0.012  sec  (post-stall)  and  0.033  sec  (pre¬ 
recovery).  However,  it  was  not  possible  to  set  the  S/RS  boundary  and  simulate  the  actual 
surge  frequency.  In  fact,  there  were  no  values  for  Xp  which  caused  the  simulation  to 
exhibit  die  actual  surge  frequency  at  100%  discharge  volume  and  81%  compressor  peed. 
The  possible  significance  of  this  will  be  discussed  later.  Once  the  values  for  Xp  were 
found,  they  were  kept  constant  in  all  simulations  depite  changes  to  the  other  variables. 

Discharge  plenum  volume  in  the  model  was  varied  by  increasing  or  decreasing  the 
size  of  several  elemental  control  volumes  in  the  rearward  section.  Since  the  elemental 
control  volumes  were  sometimes  altered  a  great  deal,  care  was  taken  so  that  the  stability 
criterion  was  not  violated  by  any  of  the  volumes.  The  model  time-stp,  At,  was  to  be 
kept  constant  in  all  cases,  so  the  time  required  for  a  sound  wave  to  travel  from  one  side 
of  a  control  volume  to  another  had  to  be  greater  than  or  equal  to  At.  This  meant  that  the 
length  of  any  particular  elemental  control  volume  could  not  be  less  than  some  set  amount. 

The  behavior  simulated  by  the  model  using  80%,  100%  and  120%  nominal  discharge 
volume  at  79%  and  81%  conpressor  speed  was  compared  to  actual  test  data.  These  six 
cases  are  shown  in  Table  6.3.  In  all  but  one  case,  the  type  of  instability  the  model 
simulated  matched  the  type  actually  observed.  However,  at  120%  volume  and  79% 
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^)eed,  the  model  Simula^  surge  while  rotating  stall  occurred  in  actuality.  More 
in^rtantly,  tlK  surge  frequencies  from  the  simulations  were  much  higher  than  those  frmn 
the  actual  test  data.  At  diis  point,  it  was  clear  that  the  model  failed  to  simulate  actual 
conq»essor  behavior  t»ing  actual  geometry  and  operating  ctmditions. 
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Next,  the  model  was  ran  with  increasingly  larger  discharge  volumes  at  a  constant 
ctmipiessor  speed  of  81%.  As  expected  the  surge  frequency  decreased  as  the  volume  was 
increased.  Figure  6.13  ^ows  diis  trend,  where  surge  frequency  is  {dotted  against 
discharge  plenum  volume.  Near  160%  nominal  discharge  volume  die  surge  frequencies 
from  the  simulation  nearly  matched  the  surge  frequencies  observed  from  the  test  data. 
Thus,  the  model  produced  surge  frequencies  near  actual  frequencies  only  by  using  larger 
volumes  than  those  presumably  tested.  Hence,  it  seemed  possible  that  the  discharge 
volumes  used  during  the  test  wae  actually  larger  than  designed. 

In  all  simulations  discussed  so  far,  die  coit^iressor  time  constants  were  held  constant. 
Therefore,  when  the  model  was  run  using  160%  volume  at  81%  qpeed,  the  simulation  did 
not  match  actual  compressor  behavior,  even  though  the  frequencies  were  neady  the  same. 
This  was  because  the  simulation  did  not  stall  at  79%  ^leed  as  the  compressor  did  in 
reality.  In  fact,  when  using  160%  simulated  volunre  the  model  did  not  exhibit  rotating 
stall  unless  the  simulated  conqitessor  speed  was  considerably  below  79%.  This  was  the 
result  of  using  the  same  Xp  values  in  all  cases. 

In  order  to  simulate  actual  compressor  behavior  using  larger  volumes,  it  was 
necessary  to  recalibrate  the  model.  This  was  done  by  a  trial-and-error  process  where 
different  volume  and  Xp  values  were  used.  After  picking  stnne  values  for  the  volume, 
values  for  Xp  were  found  which  caused  the  simulation  to  surge  at  81%  speed  and  rotating 
stall  at  79%  q>eed.  Thai  the  simulation  surge  frequency  was  conq}ared  to  the  actual 
surge  frequency  at  100%  nominal  volume.  Eventually,  a  simulation  volume  was  found 
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Table  6.3  Model  and  actual  data  comparison 
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where  both  the  surge  finquency  and  S/RS  boundary  matched  actual  test  data. 

Using  190%  simulated  volume,  the  model  exhibited  simulated  surge  behavior  very 
similar  to  that  actually  observed  at  100%  test  volume.  The  corresponding  values  of  Xj, 
were  0.012  seconds  (post-stall)  and  0.017  seconds  (pre-recovery),  hr  Figure  6.14, 
simulation  results  are  compattd  to  actual  test  data.  Both  pressure  transients  shown  are 
from  the  con^ressor  discharge.  As  the  figure  shows,  the  simulation  and  actual  pressure 
transients  neady  overlap.  If  the  model  is  a  correct  representation  of  the  compressor,  the 
simulation  results  seem  to  suggest  that  the  compressor  was  affected  by  more  volume  than 
was  present  in  die  discharge  plenum. 

Assuming  that  there  was  more  volume  than  designed,  the  compressor  would  then 
have  been  expected  to  exhibit  the  same  surge  ffequency  as  the  simulation.  Moreover,  the 
changes  in  volume  would  have  been  smaller  (by  percentage)  than  desired.  Instead  of 
±20%  nominal,  the  volume  changes  would  have  been  nearly  ±10%,  assuming  that  the 
actual  nominal  volume  was  in  fact  190%  of  the  design  value.  Smaller  changes  in  volume 
would  have  produced  smaller  changes  in  the  S/RS  boundary,  causing  the  actual  boundary 
curve  to  appear  "flatter."  Indeed,  the  S/RS  boundary  curve  based  on  the  actual  data 
appeared  to  be  a  flat  line,  which  could  be  explained  if  the  higher  volume  hypothesis  was 
correct.  There  is  also  the  possibility  that  some  flow  {dienomenon  occurred  which  caused 
the  plenum  volume  to  appear  larger  to  the  compressor  than  it  really  was  In  fact,  this  was 
probably  the  case,  since  an  error  as  large  as  190%  in  die  design  of  the  plenum  volume 
is  very  unlikely.  In  any  event,  if  there  was  more  volume  "in  play"  than  expected,  then 
it  had  to  come  from  somewhere.  Therefore,  to  prove  or  disprove  the  higher  volume 
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Actual  Dm:  81%  N.,  100%  Discharge  Pleauin  Vohioie 
OYNTECC  Simuladoo:  81%  N,.  190%  Discharge  Pleoum  Volume 
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igure  6.14  Dyntecc  simulation  and  actual  data  comparison 


hypothesis,  several  possible  sources  of  actual  and  effective  extra  volume  were  considered. 

6.2.3  Unchoked  Flow  through  the  Discharge  Valve 

During  the  c(Miq>tessor  test  and  in  both  models  previously  discussed,  the  discharge 
volume  was  defined  as  all  die  flow  volume  between  the  exit  of  the  last  stator  row  and  the 
inlet  of  the  discharge  valve.  Limiting  the  volume  size  to  the  discharge  valve  was  justified 
so  long  as  it  remained  choked.  Both  models  assumed  a  choked  throttle  (discharge  valve), 
and  in  most  cases  this  is  true.  However,  the  large  pressure  drop  associated  with  surge  in 
lO-stage  test  compressor  may  have  caused  brief  periods  of  unchoked  flow  in  the  discharge 
valve.  During  those  periods,  the  total  discharge  volume  may  have  included  some  of  the 
volume  behind  die  discharge  valve.  This  may  have  caused  the  surge  frequency  to 
decrease  or  even  have  caused  surge  to  occur  where  rotating  stall  was  expected. 

There  are  several  items  which  support  the  possibility  of  unchoked  flow.  The  volume 
downstream  of  the  discharge  valve,  which  was  made  up  of  the  discharge  ducting  and 
mass  flow  measurement  venturi,  had  a  Helmholtz  (or  organ  pipe)  frequency  of  S  to  10 
Hz.  This  was  neariy  the  same  as  the  observed  surge  frequency,  which  suggests  that  the 
discharge  duct  volume  was  affecting  surge  behavior.  Furthermore,  the  minimum 
pressures  during  surge  blowdowns  were  on  tlw  order  of  130  kPa,  nearly  atmospheric. 
Such  low  pressures  could  cause  the  discharge  valve  to  become  unchoked,  or  even  cause 
reverse  flow  if  the  discharge  ducting  does  not  depressurize  as  fast  as  the  plenum. 
Evidence  of  unchoked  flow  through  the  discharge  valve  exists  in  surge  data  from  the  E’ 
compressor  test  [7],  so  the  phoiomenon  is  not  widiout  precedent.  Certainly,  unchoking 
the  discharge  valve  could  provide  the  excess  volume  which  the  compressor  model 


84 


simulations  suggest  existed. 

In  order  for  the  valve  to  unchoke,  the  pressure  in  the  discharge  plenum  {P^)  must 

reach  some  minimum  value,  such  that  the  relation: 

P. 

2  0J283 

is  valitf .  For  instance,  if  the  static  pressure  in  the  discharge  ducting  was  atmosi^eric 
(-  101.3  kPa),  then  the  when  the  pressure  in  the  discharge  plenum  dropped  below  -192 
kPa,  the  mass  flow  through  the  valve  would  have  become  unchoked.  Since  there  was  no 
pressure  instrumentation  immediately  downstream  of  the  discharge  valve,  there  is  no  way 
to  know  exactly  what  the  pressure  in  the  duct,  ,  really  was.  However,  it  seems  certain 
that  when  the  effects  of  friction  and  flow  constriction  at  the  venturi  are  considered,  the 
pressure  near  the  discharge  valve  must  be  higher  than  the  pressure  at  the  exit  of  the 
venturi,  which  was  nearly  atmospheric.  Referring  to  the  discharge  plenum  pressure  data 
in  Appendix  C,  near  atmospheric  pressure  in  the  discharge  ducting  would  have  caused  the 
valve  to  become  unchoked  for  brief  periods  during  blowdown  and  repressurization. 
Discharge  duct  pressures  higher  than  atmospheric  would  have  caused  the  valve  to  become 
unchoked  for  even  longer  periods. 

An  attempt  was  made  to  estimate  the  pressure  difference  between  the  venturi  inlet  and 
the  discharge  valve  exit  by  determining  the  pressure  loses  due  to  friction  in  the  discharge 
duct.  Assuming  the  duct  wails  were  smooth  and  using  the  nominal  dimensions  of  the 
duct  system,  the  effects  of  friction  were  found  to  be  nearly  negligible  due  to  the  relatively 
low  flow  velocities  (Mach#  <  0.1)  near  stall.  For  this  reason,  the  pressures  near  the 

’  For  a  shaip-adged  orifice,  the  choking  pressure  ratio  may  be  greater. 
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discharge  valve  exit  were  assumed  to  be  nearly  equal  to  those  at  die  venturi  inlet. 

Time-averaged  data  taken  near  stall  at  81%  speed  indicated  that  the  pressure  and  mass 
flow  at  the  venturi  inlet  were  ~110  kPa  and  -10  kg/sec,  respectively.  However,  high- 
lesponse  data  obtained  befine  and  during  surge  indicated  that  the  discharge  valve  area, 
and  hence  the  mass  flow,  was  less  than  the  near-stall  time-averaged  value.  Because  of 
this  discrepancy,  an  effort  was  made  to  estimate  the  duct  pressure  and  discharge  valve 
Mach  number  at  stall  inception  based  on  the,  high-response  data.  A  simple  one¬ 
dimensional  model  based  on  continuity  principles  and  assuming  ismtropic  flow  conditions 
was  used  to  estimate  die  duct  pressure  and  valve  Mach  number  based  on  discharge 
plenum  data.  Since  the  flow  velocities  in  the  duct  were  low  except  in  the  discharge 
valve,  total  and  static  pressures  were  assumed  to  be  qiproximately  equal  in  all  cases. 
This  model  was  coded  in  forikan  for  use  on  a  computer. 

High-response  pre-stall  pressure  data  were  used  as  model  irqiut  data  in  order  to 
estimate  the  steady-state  flow  conditions  in  the  discharge  duct  just  before  surge.  Further, 
a  sinusoidal  pressure  variation  r^roximating  surge,  shown  in  Figure  6.15,  was  used  to 
provide  input  conditions  to  the  model  in  order  estimate  the  transient  flow  conditions  in 
the  discharge  duct  during  surge.  The  resulting  duct  pressure  and  Mach  number  estimates 
are  shown  in  Figure  6.16. 

The  qiproximate  amount  of  time  the  discharge  valve  was  unchoked  during  surge  was 
determined  by  assuming  diat  the  model  estimates  for  the  pre-stall  and  average  post-stall 
duct  pressures  were  correct.  According  to  the  model,  duct  pressure  varied  between  104 
and  107  kPa.  Based  on  this  pressure  range  and  the  ideal  choking  pressure  ratio,  the 
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discharge  valve  became  unchoked  whenever  the  discharge  plenum  pressure  dropped  below 
between  197  to  203  kPa  during  surge.  Pressure  data  firom  the  seven  surge  events  that 
occurred  at  81%  q)eed  were  conq>ared  to  this  minimum  pressure  range.  The  average 
length  of  time  that  the  valve  was  unchoked  based  on  all  seven  surge  evoits  was  23 
milliseconds  per  cycle,  which  was  an  average  of  11.1  percent  of  the  surge  cycle  period. 

There  was  also  a  variation  in  the  length  of  time  the  valve  remained  unchoked 
dqrending  on  the  plenum  volume.  At  80%,  100%,  120%  plenum  volumes,  the  unchoked 
time  lengths  were  21,  23,  and  28  milliseconds,  respectively.  This  is  a  reflection  of  how 
long  it  took  the  compressor  to  repressurize  the  plenum  in  each  case,  and  more 
in^rtantly,  demonstrates  that  the  plenum  volume  had  at  least  a  small  effect  on  the  surge 
behavior. 

Whether  or  not  the  estimated  23  milliseconds  of  unchoked  flow  during  each  cycle  had 
any  effect  on  the  instability  behavior  of  the  compressor  remains  an  unanswered  question. 
Even  if  some  of  the  discharge  duct  volume  did  contribute  to  the  overall  discharge  volume, 
it  is  difficult  to  determine  exactly  how  much.  Without  a  region  of  choked  flow,  defining 
the  boundaries  of  the  volume  to  be  included  in  any  sort  of  surge  behavior  estimate 
becomes  arbitrary.  It  is  also  important  to  remember  that  the  calculations  made  above  are 
based  on  several  assunqrtions  which  may  or  may  not  be  correct.  For  instance,  the  static 
pressure  in  the  exit  flow  of  the  discharge  valve  was  assumed  to  be  nearly  equal  to  the 
static  pressure  immediately  behind  the  valve.  In  reality,  the  static  pressure  in  the  flow 
is  probably  lower,  which  would  cause  the  valve  to  remain  choked  for  a  longer  period  of 
time.  However,  it  was  possible  to  form  an  hypothesis  concerning  the  S/RS  boundary 
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Figure  6.16 


assuming  Aat  the  valve  does  become  unchoked  for  a  significant  length  of  time. 
Unchoked  discharge  valve  hypothesis 

When  the  conqnessor  wait  into  eitha  surge  or  rotating  stall,  the  initial  behavior  of 
both  i^ienomena  were  identical.  In  bodi  cases,  the  plenum  ^pressurized  nqiidly  as  a 
result  of  conqnessor  inability  to  maintain  a  pressure  rise  across  itself.  Surge  and  rotating 
stall  behavior  were  nearly  indistinguishable  until  the  plenum  began  to  repressurize.  At 
this  point,  if  the  conqiressor  recovered  and  began  to  operate  on  the  unstalled  pressure 
characteristic  (see  Figure  1.4),  anotha  surge  cycle  began.  If  the  conqiressor  was  unable 
to  recova  (Figure  1.3),  the  compressor  appeared  to  "almost  surge"  but  then  failed  to  reach 
peak  pressure  and  slipped  into  rotating  stall. 

Since  the  behavior  of  the  two  phenomena  are  similar  until  repressurization,  it  is 
reasonable  to  assume  that  there  was  no  difference  between  surge  and  rotating  stall  until 
the  conqnessor  began  to  repressurize.  The  amount  of  volume  in  the  discharge  plenum 
was,  therefore,  irrelevant  at  first.  Assuming  that  the  characta  of  the  instability  (surge  or 
rotating  stall)  was  daeimined  just  as  the  compressor  began  to  repressurize,  the  total 
discharge  volume  became  important  only  at  this  point. 

As  the  compressor  tried  to  recover,  it  began  to  repressurize  the  volume  behind  it.  If 
the  discharge  valve  was  choked,  the  compressor  communicated  only  with  the  volume  up 
to  the  discharge  valve,  which  was  the  plenum  volume.  Howeva,  if  the  valve  was 
unchoked,  the  conq>ressor  communicated  with  the  volume  in  the  plenum  and  in  the 
discharge  ducting.  The  total  volume  including  the  ducting  was,  of  course,  much  larger 
dian  the  plenum  volume  alone.  If  the  valve  was  unchoked,  ±20%  changes  in  plenum 
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volume  w<nild  have  been  almost  insignificant  compared  to  the  total  volume.  In  fact,  these 
comparatively  small  changes  in  volume  would  have  been  nearly  unnoticed  by  the 
c(»npressor.  The  effect:  little  apparent  variation  in  instability  behavior  based  on  plenum 
volume  changes. 

Once  the  conpressor  began  to  surge  or  stall  based  (m  the  c(»nbined  volume,  it  was 
"locked"  into  that  behavior.  This  was  so  because  when  the  discharge  valve  again  became 
choked,  the  compressor  had  already  either  recovered  or  not  recovered.  Therefore,  the 
plenum  volume  alone  was  not  a  factor  in  determining  whether  the  compressor  surged  or 
went  into  rotating  stall.  However,  the  plenum  volume  did  help  determine  the  surge 
frequency,  since  the  compressor  had  to  repressurize  only  the  plenum  volume  once  the 
valve  again  became  choked. 

According  to  a  lumped-parameter  model  run  by  Ganuu;he  [21],  a  rotating  stall  event 
can  be  initiated  during  die  near  zero  flow  portion  of  a  surge  cycle  by  means  of  a  forcing 
function.  Zero  mass  flow  occurs  just  before  the  minimum  blowdown  pressure,  which  is 
almost  precisely  where  the  hypothesis  says  the  delineation  between  surge  and  rotating 
stall  was  made.  So  pedups  it  is  trae  that  the  character  of  a  conpressor  instability  event, 
whether  it  is  surge  or  rotating  stall,  is  definable  at  a  minimum  of  two  locations  on  the 
conpressor  characteristic:  near  stall  and  near  zero  flow. 

It  remains  to  be  seen  if  this  hypothesis  can  be  proved.  Since  there  was  no  pressure 
data  available  immediately  behind  the  discharge  valve,  it  was  impossible  to  say  for  certain 
that  the  discharge  valve  was  indeed  unchoked  for  the  estimate  period  of  time.  Extremely 
brief  periods  of  unchoked  flow  would  probably  not  have  affected  conpressor  instability 
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bdiavior,  so  it  is  critical  to  tfiis  hypothesis  that  the  valve  be  unchoked 
for  smnc  unspecified  minimum  amount  of  time.  Further,  there  was  no  way  to  be  sure  at 
what  point  compressor  instability  behavior  was  irrevocably  defined.  Therefore,  this 
l^pothesis  remains  interesting  but  lume-the-less  uiqnoven. 

6.2.4  Combustor  insert  filler  volume 

In  the  actual  engine,  the  volume  immediately  behind  the  compressor  discharge  is 
taken  up  by  the  combustor.  Unfortunately,  it  was  difficult  to  design  variable  volume 
hardware  which  would  allow  the  combustor  volume  to  vary  without  upsetting  the  flow 
path.  To  solve  this  problem,  the  test  engineers  made  the  assurrption  that  flow  velocity 
inside  the  combustor  itself  was  low,  which  allowed  them  to  replace  the  volume  inside  the 
combustor  with  an  open  plenum  volume  just  aft  of  the  combustor  location.  At  low  flow 
velocities,  it  would  not  matter  where  the  volume  was  located.  Therefore,  the  variable 
volume  hardware  could  be  placed  in  this  alternate  location  widiout  upsetting  the  flow 
path. 

The  combustor  flow  path  was  maintained  by  the  incorporation  of  a  hollow  insert 
(filler  volume)  as  shown  in  Figure  3.4.  Four  6.3S-mm  holes  were  drilled  through  the 
insert  walls  to  allow  air  to  circulate  between  the  filler  volume  and  the  discharge  plenum 
volume.  These  holes  were  included  so  that  during  steady-state  operation  the  pressure 
inside  the  insert  could  equalize  widi  the  pressure  in  the  discharge  plenum.  This  reduced 
the  risk  of  pressure  induced  stresses  which  could  cause  deformation  or  damage. 

Although  interaction  between  the  volumes  was  desirable  during  steady-state  operation, 
it  was  not  during  transient  operation.  The  small,  6.3S-mm  hole  size  was  originally  chosen 
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so  dut  during  transient  operation,  the  volume  inside  the  insert  wmild  not  interact  widi  the 
discharge  plenum  volume.  To  verify  this,  a  simple  flow  model  was  written  to  see  if  die 
air  mass  inside  the  combustor  could  have  significandy  cmtributed  to  the  overall 
compressor  mass  flow  during  a  pressure  transient.  The  model,  which  was  coded  in 
FORIItAN  for  computer  use,  was  fundamentally  the  same  model  used  to  estimate  flow 
behavior  in  die  discharge  ducting. 

The  same  sinusoidal  surge  representation  and  initial  conditions  used  in  the  discharge 
duct  model  were  used  for  the  iQller  volume  model.  Unlike  the  duct  volume  model,  where 
the  direction  and  amount  of  mass  flow  through  the  flow  tunings  were  not  necessarily 
the  same  or  equal,  the  geometry  of  the  insert  guaranteed  that  the  mass  flows  through  the 
four  holes  would  be  equal  and  in  the  same  direction  at  all  times.  Thus  the  flow  areas  of 
all  four  holes  were  summed  so  that  the  model  could  treat  all  four  holes  as  a  single 
opening.  For  the  purposes  of  comparison,  the  model  was  also  run  using  a.  step  drop  in 
pressure  as  an  itq)ut.  This  was  done  to  (tetermine  how  long  the  filler  volume  would  take 
to  blow  down  with  a  constant  imposed  plenum  pressure  as  0|^sed  to  a  varying  one.  For 
consistency,  the  minimum  pressure  in  the  step  pressure  drop  was  set  equal  to  the  average 
pressure  of  the  sinusoidal  surge  representation. 

The  results  from  the  model  are  shown  in  Figures  6.17  and  6.18.  In  Figure  6.17,  filler 
volume  pressure  and  mass  flow  are  shown  as  a  function  of  time.  As  expected,  the 
pressure  fell  exponentially  toward  some  minimum  value  and  the  mean  mass  flow 
iq)proached  zero  during  die  simulated  surge.  Even  at  its  peak  value,  the  mean  filler 
volume  mass  flow  was  far  smaller  than  the  total  mass  flow.  According  to  the  model,  the 
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mean  mass  flow  out  of  the  filler  volume  never  exceeded  0.2%  of  the  total  mass  flow. 
This  is  shown  in  Figure  6.18.  The  blowdown  period  observed  during  simulated  surge  was 
also  long  (  >15  seconds)  when  con^axed  to  that  observed  during  the  step  drop  in 
pressure.  The  typical  time  length  of  a  recorded  transient  from  this  test  was  less  than  2 
seconds,  so  little  air  was  likely  to  have  escaped  from  the  filler  volume  during  that  time. 
It  is  unlikely  that  this  small  amount  of  air  contributed  significantly  to  the  instability 
behavior  of  the  test  conq)ressor,  so  it  was  possible  to  coiiclude  that  the  filler  volume  had 
little  or  no  effect  on  the  S/RS  boundary  in  this  conq>ressor  based  on  these  results. 

6.2.5  Other  possible  sources  of  extra  discharge  volume 

The  possibility  of  other  sources  of  volume  in  and  around  the  compressor  discharge 
was  investigated  by  closely  examining  the  design  drawings  and  querying  those  personnel 
who  actually  worked  with  the  test  compressor.  Neither  the  drawings  nor  the  experiences 
of  those  questioned  could  account  for  any  significant  amount  of  extra  discharge  volume. 
A  certain  amount  of  volume  related  to  the  seals  and  coiuiections  to  the  discharge  ducting 
existed,  but  these  were  clearly  negligible.  There  were  no  other  sources  of  extra  discharge 
volume  other  than  the  filler  volume  and  discharge  ducting  as  already  discussed  in 
Sections  6.2.4  and  6.2.5. 

6.2.6  Choked  flow  In  the  rear  stages  of  the  compressor 

Earlier  in  this  report,  the  stage-by-stage  compressor  model  DYNIECC  was  used  to 
demonstrate  the  possibility  of  extra  volume  afrecting  the  S/RS  boundary  behavior  of  this 
compressor.  One  of  the  basic  assumptions  which  allowed  this  model  to  be  used  was  that 
iMither  sonic  nor  supersonic  flow  occurred  anywhere  in  the  model  except  at  the  discharge. 
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volume  dynamic  behavior  estimate 
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Figure  6.18  Insert  volume  mass  flow  ratio 


In  fact,  if  the  iiqmt  getmietry  were  such  that  sonic  flow  occurred,  the  model  would  have 
failed  utteily. 

In  his  technical  report,  Copenhaver  [2]  suggested  that  tlM  unblocked  portitm  of  the 
temh  stage  rotor  was  probably  choked  during  rotating  stall.  In  the  model,  this  is 
accounted  for  in  the  tenth  stage  characteristics.  When  the  surge  data  used  in  this  report 
were  reviewed,  the  possibility  of  choked  flow  in  the  exit  stator  row  was  also  discovered. 
E>ue  to  the  stator-rotor  stage  definition  used  during  the  compilation  of  the  stage 
characteristics,  the  flow  conditions  in  the  exit  stator  row  may  not  have  been  included  in 
the  tenth  stage  characteristics.  Furthermore,  dyniecx:  does  not  model  flow  through  blade 
rows,  so  it  is  possible  that  instances  choked  flow  occurred  in  the  exit  stator  row  of  the 
actual  conq>ressor  that  were  not  accounted  for  in  the  model. 

Even  if  the  tenth  stage  characteristics  included  exit  stator  flow  conditions,  another 
assumption  made  in  order  to  use  the  DYNTECC  model  could  also  be  invalidated  by  the 
possibility  of  choked  exit  stator  flow.  When  compiling  the  stage  characteristics  for  use 
with  the  DYNTECC  model,  it  was  necessary  to  assume  that  the  stall  characteristics  for 
78.5%  speed  could  be  substituted  for  those  at  81%  speed.  Since  die  compressor  surged 
at  81%  speed,  and  the  researchers  responsible  for  the  compressor  did  not  want  to  risk 
causing  damage  to  the  compressor  or  the  facility,  steady-state  data  were  not  obtained 
under  those  operating  conditions.  If  the  exit  stators  did  choke  during  surge  but  not  during 
rotating  stall,  this  could  have  caused  the  78.5%  and  81%  speed  stall  characteristics  to 
differ  greatly.  This  would  invalidate  the  assumption  of  similar  stall  characteristics  at 
78.5%  and  81%  speed. 
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To  explore  the  possibility  of  choked  flow  in  the  stators  during  surge,  a  simple 
isentropic  flow  model  was  devised.  Using  the  high-response  static  and  total  pressure  data 
from  die  conqnessor  discharge,  the  Mach  number  in  the  exit  stator  row  was  estimated. 
The  flow  area  through  the  stator  row  was  critical  in  these  estimates.  Since  die  flow  area 
(tepends  entirely  upon  the  flow  angle  through  d^  exit  stator  row,  three  different  values 
for  the  flow  area  were  run.  These  three  different  values  represented  the  effect  three 
different  flow  angles  had  on  the  flow  area  in  the  stator  row.  Boundary  layer  effects  were 
ignored  in  these  calculations,  so  the  effective  areas  are  smaller  than  those  used  in  this 
model.  This  would  cause  the  stators  to  become  choked  more  often  than  the  model 
predicted.  Some  of  the  model  results  are  shown  in  Figure  6.19  and  Table  6.4. 

According  to  the  model,  flow  through  the  stators  may  have  been  choked  for  time 
periods  of  up  to  16  milliseconds.  These  possible  periods  of  choked  flow  always  occurred 
during  the  repressurization  phase  of  a  surge  cycle.  Shorter  time  periods  were  expected 
in  those  cases  where  the  plenum  volume  was  smaller,  since  it  takes  less  time  to 
repressuxize  smaller  volumes.  The  apparent  disagreement  with  this  trend  at  120%  volume 
is  probably  not  significant;  only  three  surge  cycles  at  that  speed  and  volume  were 
available  for  averaging.  Although  these  data  are  interesting  in  a  qualitative  sense,  it  is 
important  not  to  attach  too  much  significance  to  these  calculations.  Unfortunately,  there 
was  insufficient  data  available  to  determine  the  flow  angle  in  the  exit  stator  row,  so  there 
was  no  way  to  prove  that  the  stators  were  choked  or  unchoked  during  any  of  the 
documented  surge  or  rotating  stall  events. 


98 


6.2.7  Other  possible  reasons  for  S/RS  boundary  behavior 

Several  other  possible  reastxis  for  the  unexpected  S/RS  boundary  behavior  abound. 
One  item  that  could  have  had  a  significant  affect  on  the  instability  behavior  of  this 
conqiressor,  but  was  almost  con^letely  ignored  in  diis  report,  was  the  tq>stream  volume. 
The  DYNlECC  noodel  included  upstream  volume  up  to  the  bellmouth,  but  there  was  a  large 
amount  of  volume  in  front  of  the  bellmouth,  including  the  flow  conditioning  barrel,  that 
went  ignored.  In  the  case  of  the  B  parameter  calculations,  virtually  none  of  the  volume 
in  front  of  the  conqnessor  was  considered.  Only  volume  up  to  -0.20  meter  in  front  of 
the  inlet  guide  vanes  was  included  in  the  form  of  the  effective  conqnessor  length, 
Exactly  how  much  of  the  upstream  volume  should  have  been  used  in  the  instability 
behavior  estimates  was  not  clear.  It  was  decided  that  much  of  the  forward  volume  could 
be  ignored  since  the  flow  velocities  upstream  of  the  beUmouth  were  relatively  low.  It  is 
the  feeling  of  this  researcher  that  omitting  much  of  the  upstream  volume  was  justified, 
but  it  is  certainly  possible  that  much  more  of  the  upstream  volume  should  have  been 
considered  in  both  the  stage-by-stage  model  and  in  the  B  parameter  calculations. 
Determining  how  much  of  this  volume  to  use  is  a  recommendation  for  future  research. 

As  previously  mentioned,  the  simple,  lun^d-parameter  model  used  to  estimate  the 
S/RS  boundary  does  not  account  for  die  compressibility  of  air  across  the  compressor 
itself.  It  is,  therefore,  possible  that  had  this  been  somehow  included,  the  estimated 
behavior  would  have  matched  the  observed  behavior.  This  is  doubtful,  since  the  S/RS 
boundary  didn’t  seem  to  move  at  all.  It  seems  likely  that  even  including  compressibility 
the  lumped-parameter  model  would  have  predicted  some  shift  in  the  SA^S  boundary, 
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although  this  is  not  necessarily  so. 

After  conq}leting  this  research  effort,  it  was  parent  that  the  S/RS  boundary 
bdiavior  could  not  be  predicted  with  existing  theories  unless  certain  assun^ons 
axKeming  die  flow  volume  were  made.  Therefore,  it  seems  likely  diat  assun^ions 
ctmceming  the  flow  volume  are  the  reason  for  the  disagreement  between  observed  and 
predicted  S/RS  boundary  behavior. 


7.0  SUMMARY  AND  CONCLUSIONS 

This  report  presented  some  of  die  surge  data  associated  widi  a  10-stage  core 
ccxi^iressor  test  cmiducted  at  die  CRF.  Of  the  at  least  24  documented  instances  of  surge 
observed  during  this  test,  12  were  used  in  this  report.  Pressure  data  from  three  difrerent 
axial  locations  in  the  compressor  were  used  to  help  characterize  the  surge  events.  It  was 
clear  from  this  data  that  the  surge/rotating  stall  bourulary  did  not  move  significandy  whoi 
the  discharge  plenum  volume  was  changed  ±20%.  Possible  reasons  for  this  lack  of 
movement  were  discussed,  and  the  results  suggest  that  care  must  be  taken  when  qiplying 
current  theories  to  the  CUF  conpvssor. 

The  following  are  observations  and  conclusions  concerning  the  general  character  of 
the  12  surge  events  used  in  this  research  effort: 

1.  The  12  surge  events  documented  in  this  report  are  exanqiles  of  classic  surge.  The 
ri^id  blowdown  and  slower  repressurization  periods  observed  in  each  surge  cycle  are 
typical  of  this  type  of  surge.  Further,  the  surge  frequencies  of  -S  to  -7  cps  are  small 
conqiared  ^  >  the  rotating  stall  frequencies  of  -50  to  -100  cps,  which  were  also  typical. 

2.  Periods  of  reverse  flow  during  the  blowdown  period  of  each  surge  cycle  were 
evident.  Data  obtained  from  pressure  instrumentation  in  the  conqjtessor  discharge 
demonstrate  this  conclusively. 

3'  Plots  of  discharge  pressure  as  a  function  of  time  showed  that  rotating  stall  occurred 
simultaneously  with  surge  during  die  blowdown  periods  of  all  12  surge  events.  If  they 
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had  ctmtinued  past  the  blowdown  period  to  become  fully  devel(^)ed,  the  stall  cells  would 
have  exhibited  passage  frequencies  of  70  to  100  cps,  which  aie  ccmqdetely  leastmable  for 
the  speeds  at  whidi  they  occuned. 

4.  The  compressor  was  tested  using  ditee  different  discharge  fdenum  vdiume  values: 
80%,  100%  and  120%  nominal  volume.  The  surge  firequoicy  may  have  changed  slightly 
when  the  plenum  volume  was  changed.  There  is  insufficient  data  to  prove  this,  and  evm 
if  such  a  change  occurred,  it  was  very  slight  indeed  (from  -4.3  to  -4.8  cps  at  81%  speed). 

The  following  are  observations  and  theories  concerning  the  bduivior  of  foe 
surge/rotating  stall  boundary: 

1.  Theory  predicts  a  shift  in  tlw  surge/rotating  stall  boundary  when  the  plenum 
volume  is  altered,  but  foe  test  con^ressor  showed  no  such  shift.  In^ed,  the 
surgeAotating  stall  speed  boundary  seemed  to  remain  constant  de^ite  ploium  volume 
changes. 

2.  Six  theories  which  attempt  to  explain  foe  constant  surge/rotating  stall  boundary 
were  formulated  and  investigated  in  this  report.  Table  7.1  lists  these  theories  and  the 
results  of  the  investigations  concerning  their  merit.  No  theory  was  proved  to  be  the  cause 
of  the  surge/rotating  stall  boundary  bdiavior.  However,  foe  unchdced  discharge  valve 
theory  was  found  to  be  the  most  reasoii^ble  explanation.  After  reviewing  foe  research 
conducted  with  other  high-pressure  test  conqrressors,  such  as  foe  E^,  it  seems  possible  that 
this  sort  of  behavior  is  not  unique  to  this  compressor.  If  true,  this  theory  and  explanation 
differentiates  foe  surge  behavior  of  low-  and  high-pte:^sure  compression  systems. 
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Table  7.1  Surge/rotating  stall  boundary  behavior  theories 


Theoiy 


UnelMfcwifkyw 

veiv«  allowed  extra 
vohnni  lobe 
perceived  by  lha 


Possible  Affect 
on  Instability 
Behavior 


Ibcwe  ia  likaHlMad  of 
larfB  Umtmd  of  ratatBg 
mJI:  rodaetioa  ia  twf» 


Mayakcrtfwi 


DYNTEOC  10  limalaH 
S/RS  boonrtfy  bahavior 


votaM  JaapiiB  baing 
daaigaad  aot  lodo  ao 


• - i_  trvainirtrtft 

ncieeaB  tn  uxeiinoeo  or 
aarga  ioMead  of  raiaUng 
atatl;  radaelioB  bi  aarga 
fragaancy 


Moaa  of  dw  bibt 
vohraaa  ahoaU  ba«o 
oeeii  ceenoBieQ  n  ew 


comribaiod  to  tha 


hcioaae  in  HkeHliood  of 
aarga  iitaari  of  rattling 
attll;  radactiaa  in  aaiga 
fiai|aancy 


Proven 

or 

Diq)coved 

Actual 

Affect 

Neitber 

(Uncbohad  flow  did 
occar.  bat  poaaably  for 
only  very  brief  parioda) 

Uoknown 

Nddier 

Unknown 

Diapioved 

None 

Neither 

(Not  raaearabed) 

Unknown 

Disproved 
(No  extra  volame  from 
aeala  or  cannectiona) 

None 

Neitber 
(Not  raaaarched) 

Unknown 
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3.  Two  of  the  theories  presented  in  Table  7.1  were  not  researdied,  but  were 
formulated  as  a  result  of  this  research  effort.  They  are  also  listed  am<mg  the 
recommendations  as  topics  for  future  research. 
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8.0  RECOMMENDATIONS 


Much  information  ccmceming  the  surge  behavior  of  a  high-q)eed,  multi-stage 
conquessor  has  been  presented  in  diis  report.  Mote  such  information  should  be  generated 
to  help  continue  the  investigation  of  both  the  surge/rotating  stall  boundary  behavior  of  this 
con^ressor  and  of  the  surge  behavior  of  high-speed  compressors  in  general.  Certainly, 
a  follow-on  test  of  this  or  a  similar  con^ressor  is  recommeixled,  with  some  of  the 
following  in^rovements  included. 

1.  High-response  pressure  and  temperature  instrumentation  should  be  placed  at  the 
entrance  and  exit  of  the  discharge  valve.  Similar  instrumentation  should  be  placed  at  the 
entrance  and  exit  of  both  the  exit  stator  row  and  the  tenth  stage.  This  would  allow 
researchers  to  determine  if  choking  is  occurring  at  the  two  locations. 

2.  Surges  of  Itmger  duration  should  be  induced  and  recorded.  Further,  the  discharge 
valve  area  should  be  held  constant  for  a  significant  portion  of  those  surge  events.  The 
fundamental  problem  with  doing  this  is  that  stresses  in  the  rotor  may  become  too  extreme 
if  die  compressor  remains  in  surge  for  too  Itnig.  However,  a  longer  surge  event  with 
constant  discharge  valve  area  would  allow  a  more  confident  estimate  of  the  surge 
frequency  to  be  made.  Comparing  individual  surge  events  with  one  another  and  with 
model  generated  data  would  also  be  easier. 

3.  More  steady-state,  high-reqxmse  data  should  be  recorded  in  conjunction  widi  the 
time-averaged  data  points.  During  this  test,  researchers  sometimes  turned  off  the  analog 
tape  recorders  during  steady-state  data  points  in  an  effort  to  ctmserve  tape.  Unfortunately, 
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this  made  calibratitm  of  the  high-re^xmse  in^nimentation  more  difficult.  Without  the 
time<averaged  data,  calibration  of  the  high-re^nse  instruments  depends  on  either  the 
bench  calibratkMis  or  the  individuai  monitor  points,  both  of  whidi  are  less  accurate  tluoi 
using  the  time-averaged  data.  It  was  fortunate  that  diis  teseardier  found  sufficient  time- 
avetaged^ugh-te^Kmse  data  ovetlrp  to  calibrate  foe  high-response  instrumoits,  but  these 
calibration  points  sometimes  numbered  mily  three  or  four  points  in  a  day.  More 
calibration  points  ate  desirable.  The  more  calibration  points  there  are,  foe  more  accurate 
the  calibration  equations  become. 

Even  if  more  surge  data  ate  not  obtained  for  this  cotrqnessor,  there  are  several  items 
which  should  receive  more  attention  if  research  concerning  the  surge/rotating  stall 
boundary  of  this  compressor  ctmtinues. 

1.  Consideration  should  be  given  to  including  mote  inlet  volume  in  both  foe  lumped- 
parameter  and  foe  stage-by-stage  models.  The  additiorud  volume  may  help  account  for 
the  surge  behavior  of  this  compressor. 

2.  The  ability  to  model  choked  flow  in  the  compressor  blade  passages  should 
somehow  be  incorporated  imo  DYNTECC. 
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APPENDIX  A:  Compressor  Research  Facility 

Description 


^)pendix  A  briefly  describes  the  layout  and  operation  of  the  CRF.  Hguie  A.1  shows 
the  buildings  that  comprise  the  CRF  while  Figure  A.2  shows  a  schematic  of  the  test 
chamber. 
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The  data  presented  in  diis  report  were  obtained  during  a  compressor  te^  carried  out 
at  the  Conopressor  Research  Facility  (CRF)  at  Wright-Patterson  AFB.  A  rough  layout  of 
this  facility  is  shown  in  Figure  A.I.  The  operations  building  ccmtained  the  test  chamber, 
signal  ccmdititming  equqnnent,  and  compressor  drive  motors.  The  electrical  power 
conditioning  building  contained  such  equqnnait  while  the  test  building  contained  the 
computer  and  test  compressor  control  systems. 

The  test  was  controlled  from  the  facility  control  room  located  in  the  test  building. 
Commands  to  carry  out  all  test  procedures  were  sent  from  the  control  room  via  a  monitor 
computer,  located  in  the  copiputer  room,  to  four  control  computers  located  in  the  signal 
conditioning  room.  The  control  computers  provided  input  to  the  electrical  power 
conditioning  equipment  which  powered  the  drive  motors.  The  control  conq>uters  also 
controlled  services  that  support  variable  geometry  devices  on  the  test  compressor. 

The  test  compressor  was  driven  by  22.37-MW  (30,000  hp)  synchronous  electric  drive 
motor.  The  test  conqttessor  speed  range  was  provided  through  two  speed  increasing  gear 
boxes.  A  speed  range  of  3,000  to  15,0(X)  RPM  was  required  for  this  test  compressor. 
Within  this  range,  the  con:q}tessor  speed  was  changed  by  varying  the  frequency  of  the 
current  supplied  to  the  drive  motor.  Frequency  variation  of  the  supplied  current  was 
acconq)lished  by  the  electrical  power  conditioning  equipment. 

The  test  conqpressor  was  mounted  inside  a  6.1-m  (20.0  ft)  diameter  test  chamber, 
shown  in  Figure  A.2.  Atmospheric  air  passed  through  frlteting  equipment,  five  inlet 
valves,  and  a  series  of  flow  conditioning  elements  before  it  was  finally  drawn  into  the  test 
compressor.  Air  was  discharged  from  the  test  compressor  into  die  atmosphere  through 
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the  throttling  valve,  die  facility  discharge  collector,  and  discharge  ducting  (see  Figure 
A.2).  Exit  flow  from  the  test  compressor  was  monitored  by  a  48.26-cm  (19  in)  throat 
diameter  venturi  positioned  in  the  discharge  ducting  approximately  30.48  m  (100  ft) 
downstream  of  the  compressor  exit. 
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Figure  A.1  Compressor  Research  Facility  layout 


TO  ATMOSPHERE 


Figure  A.2  Test  compressor  installation  in  test  chanter 


APPENDIX  B:  Data  Acquisition 


^)peiidn  B  d^ails  the  high-iesp<Hise  and  time-averaged  data  acquisition  methods 
used  to  gather  data  during  die  conqiressor  test.  Measurement  uncertainty  is  also 
discussed.  Figures  B.l  and  B.2  show  schraiatics  of  the  high-response  and  time-averaged 
systems,  respectively. 
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During  the  CRF  ccmqnessor  test,  three  methods  of  pressure  data  acquisititm  woe 
oi^loyed:  high-response,  close-coupled,  and  time-averaged  data  acquisition.  These  three 
m^ods  acquired  data  with  frequency  responses  of  200  Hz,  70  Hz,  and  0.5  Hz 
respectively.  Since  tmly  high-response  and  time-averaged  data  are  presented  in  this 
report,  close-coupled  data  acquisition  will  not  be  covered  in  detail  here.  For  details  on 
the  closed-coupled  data  acquisition  system,  see  References  2  and  3. 

High-Response  Data 

IDuring  the  con:q)ressor  test,  52  channels  were  designated  as  high-re^nse  pressure 
measurement  channels.  Each  channel  had  a  nmninal  frequency  re^Kuise  of  200  Hz, 
which  was  adequate  to  document  the  fluctuating  pressures  expected.  Only  pressure 
fluctuations  with  fundamental  frequencies  of  less  than  50  Hz  were  anticq>ated  to  be  of 
importance.  This  is  higher  than  the  fundamental  frequencies  documented  during  the  surge 
events  discussed  in  this  report.  Only  6  of  the  52  available  high-response  pressure 
measurement  channels  were  used  in  this  report. 

The  pressure  transducers  used  to  obtain  the  high-response  pressure  measurements 
were  Kulite  model#  XCQ-25’s  (inlet),  XCQ-150’s  (discharge)  and  XTEL-190’s  (inlet 
stage  10).  Each  of  these  transducers,  with  the  exception  of  diose  at  the  inlet,  was 
contained  in  a  special  cooling  jacket.  Cooling  was  required  to  assure  survivability  and 
stability  of  the  transducers  in  the  test  environment.  Transducer  ranges  were  based  on  the 
maximum  pressures  expected  during  testing,  which  occurred  at  87  percent  corrected  rotor 
speed,  the  maximum  conq>ressor  speed  tested.  For  the  configuration  and  placement  of 
the  high-response  pressure  transducers  used  in  this  report  see  Section  3.1. 
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TEST  CHAMBER  I  SIGNAL  I  CONTROL 

:  CONDITIONING  :  ROOM 


Figure  B.1  High-response  data  acquisition  system  schematic 


A  schenuttk  of  the  hi^-re^XMise  data  acquisition  system  is  diown  in  Figure  B.l.  In 
this  system,  Pacific  Model  #8253  signal  ccmditioneis  provided  excitation  voltages  ro  the 
pressure  transduceis  and  handled  transducer  output  signal  filtering  and  anq>lificati<m 
during  testing.  Furthamore,  the  conditioners  provided  ac/dc  co«^ling  fm  the  ventage 
signals  from  the  transducers.  Preset  cut-off  filter  frequencies  of  10  Hz,  100  Hz,  1  kHz, 
and  wkte  band  were  available  witii  tiiese  condititHiers.  The  transducer  ouqmt  signals  were 
attenuated  by  a  two  pole  Bessel  filter  with  a  -3-db  attenuatiem  at  the  cut-off  frequency 
chosen.  The  cut-off  frequency  for  this  test  was  set  at  1  kHz  to  assure  die  desired 
frequency  response  of  200  Hz  would  be  obtained.  Transducer  output  was  dc  coupled,  so 
diat  the  dc  conqxment  could  be  used  in  calibration  procedures  to  determine  actual 
jxressures  at  die  measurement  locations  (see  Section  5).  Care  was  taken  to  avoid 
saturation  of  the  anqilifiers  and  the  recording  equqiment  since  the  dc  compements  of  the 
transducer  output  signals  were  being  retained. 

The  high-re^nse  pressure  measurements  were  recorded  on  three  Datatype  model 
#2808,  28-channel  recorders  configured  for  frequmey  modulation  (FM)  recording.  Each 
pressure  channel  was  assigned  to  a  track  (chanrwl)  on  the  recorder  and  recorded  on  analog 
tiqpe.  Signals  frmn  each  channel  were  recorded  at  a  center  frequency  of  73.6  kHz. 
Frequency  deviatiem  was  adjusted  to  allow  for  a  maximum  irqmt  signal  level  of  3  volts. 
Recording  qieed  was  set  at  38.1  cm/sec  (15  in/sec)  which  allowed  for  a  data  frequency 
respcMise  of  10  kHz.  This  is  well  above  the  required  200  Hz  response.  Channels  that 
required  jrfiase  correlation  for  post-test  conqiarison  were  recorded  cm  either  all  odd  or  all 
even  track  (chatmel)  numbers.  All  high-response  data  presented  in  this  report  were 
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recorded  on  odd  track  numbers  except  for  stage  9  static  pressure  data.  Also  recorded  on 
atudog  tape  were  data  signals  indicating  discharge  valve  position.  These  signals  were 
provided  by  a  potentiometer  used  to  rrxxutor  the  discharge  valve.  Like  the  high-resptmse 
jnessore  transducer  sigiuls,  the  dc  ctmqxment  of  the  potentitmieter  signal  was  r^ained  to 
allow  for  post-test  calibration. 

Along  with  the  data  signals  mentioned  above,  three  other  signals  were  recorded  on 
the  analog  ti9)es.  One  track  on  the  recorder  was  designated  for  recording  of  IRIG  A  time 
code,  which  allowed  for  timing  resoluticm  of  0.1  millisectMid.  Another  track  was 
designated  for  recording  of  the  tipe  drive  servo  signal  which  is  used  to  set  playback 
speed.  This  is  dcme  so  that  any  speed  fluctuations  in  tape  speed  will  be  duplicated  in 
playback.  A  third  track  was  used  to  record  a  standard  wave  signal,  which  can  be  used 
to  verify  phase  alignment  when  data  recorded  on  two  different  tipes  (but  recorded  at  the 
same  time)  are  conpared. 

Effort  was  made  to  verify  that  the  pressure  transducers  remained  stable  on  a  day-to- 
day  basis.  At  selected  times  during  testing  the  test  coirpressor  was  brought  to  a  fixed 
operating  point.  At  this  "checkpoint”  the  dc  levels  of  the  high-response  pressure 
transducers  were  noted  and  cotrpared  to  dc  levels  of  the  previous  test  period.  Differences 
in  day-to-day  dc  levels  were  dmeby  noted.  Transducer  instability  due  to  dc  drift  was  not 
significant  during  this  test.  However,  the  daily  "checiqwints"  did  help  locate  pressure 
transducers  that  were  not  functioning  properly.  The  high-response  transducers  used  in  the 
test  were  fragile  and  ixn  really  suited  to  the  test  envirmiment,  so  malfimctitm  was  not 
uncommon.  Faulty  transducers  of  this  type  were  often  noted  and  replaced  where  possible. 
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A  more  detailed  discussioa  of  the  high-r«q)onse  data  acquisition  sy^em  used  during 
the  ctmqtiessor  test  is  presented  in  a  pq)er  written  by  Cq)aihav«  and  Worland  [13]. 
Time-Averaged  Data 

Four  hundred  and  one  channels  were  set  aside  for  time-averaged  meararements  during 
the  compressor  test.  Time-averaged  pressure  and  temperature  data  woe  taken  to 
characterize  the  steady  (unstaUed)  and  quasi-steady  (rotating  stall)  poformance  of  the  test 
cwnpressor.  Time  averaged  data  were  obtained  in  digital  form  during  the  test  and 
recorded  on  digital  tspe.  This  allowed  for  easier  ctmversim  to  engineering  units  during 
post^MTOcessing.  Gose-coupled  data,  which  are  iK)t  presented  in  this  report,  were  obtained 
using  the  same  instrumentati<Mi  as  die  time-averaged  data.  These  two  acquisition 
techniques  differed  only  in  their  reflective  data  acquisition  rates  and  number  of  scans 
made  per  measuranent. 

Time-averaged  data  were  obtained  by  taking  a  fiecified  number  of  individual  sanqiles 
over  a  fiecific  time  period  and  then  calculating  the  arithmetic  average.  Tlw  calculated 
average,  a  voltage,  was  then  stored  and  used  to  calculate  the  pressure  or  teiiqierature  at 
the  measurement  location.  In  this  test,  30  individual  sanqiles  were  taken  over  a  time 
period  of  195  ms.  The  time-averaged  sanqile  rate  (frequency  tefxmse)  was  dictated  by 
the  maximum  rate  at  which  analog  signals  could  be  digitized  and  by  die  tube  length 
between  the  pressure  port  and  the  transducer.  This  limited  the  frequency  refionse  of  the 
time-averaged  acquisition  system  to  about  0.5  Hz.  The  use  of  30  sanqiles  to  calculate  the 
average  was  based  on  limitations  inqiosed  by  the  statistical  mediod  employed  to  calculate 
measurement  uncertainty. 
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A  schanadc  of  the  time-averaged  data  actpiisition  syston  is  shown  in  Figure  B.2. 
The  pressure  transducers  used  by  this  system  were  Scanivalve  model  ZCXT-M’s  ^ro- 
2perate-£alil»rate).  These  were  a  "more  rugged"  type  of  transducer  than  those  used  by 
tile  hi^response  syston  and,  therefore,  more  suited  to  tiiis  test  environmenL  However, 
it  was  still  necessary  to  control  the  tenyerature  of  the  ZOC  transduces  in  order  help 
guarantee  calibration,  zero  offset,  and  sensitivity  stability.  This  was  acconqilished  by  a 
closed  loop  cooling  system  with  cooling  water  maintained  at  21.1°  C  (70°  F).  Pressure 
transduce  signals  were  conditioned  by  Preston  (model  8300-XWBRC)  amplifiers.  Each 
anqilifier  was  ciqiable  of  outputting  a  10-volt  signal,  and  each  was  equqiped  with  a  three- 
pole  Butterwortii-Thomson  filter  with  selectable  cutolf  fiequencies  of  10  Hz,  40  Hz,  120 
Hz,  400  Hz,  and  1200  Hz.  The  analog-to-digital  (A  to  D)  converte  saturated  at  S  volts, 
so  appropriate  gains  were  selected  for  each  channel  to  eisure  that  the  maximum  output 
fnxn  the  anq>lifier  was  as  close  to  5  volts  as  possible.  The  cutoff  fiequency  chosen  for 
the  te^  was  10  Hz  (for  time-averi^ed  acquisitirm).  This  fiequency  was  the  lowest 
possible  filter  setting  and  reduced  the  level  of  60  Hz  noise  in  the  transducer  signals. 
After  crmditioning,  the  transducer  sigiuds  were  sampled  by  the  A  to  D  converter  to 
provide  the  30  individual  san:q)les  necessary  for  calculation  of  the  time-averaged 
measurements.  The  time-averaged  measurements  were  then  processed  by  the  Data 
Acquisition  Conqniter  (DAC)  and  recorded  cm  digital  tape  via  the  Main  Computer  (IBM 
4341). 

Performance  mapping  of  the  test  conqnessor  was  accon:q)lished  by  using  the  time- 
averaged  data  acquisition  system  to  provide  steady-state  and  quasi-steady-state  data  points 
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Figure  B.2  Time-averaged  data  acquisition  system  schematic 


for  the  nu^.  At  each  data  point,  a  time-averaged  meaairement  from  each  relevant 
transducer  was  takoi.  The  data  obtained  at  each  steady-state  data  point  were  then  used 
to  determine  die  operating  characteristics  of  the  compressor  at  that  operating  point.  This 
infotmatixm  was  in  turn  used  to  create  a  performance  map  of  die  compressor.  Three 
minutes  waiting  time  between  each  set  of  time-averaged  measurenaents  was  the  minimum 
time  delay  allowed.  This  was  done  to  assure  steady-state  or  quasi-steady-state  equilibrium 
had  been  reached  after  each  operating  condition  change. 

In  this  report,  time-averaged  data  in  the  form  of  overall  compressor  and  individual 
stage  characteristics  was  used  in  the  mathematical  models.  Hme-averaged  data  were  also 
used  to  convert  high-response  transducer  data  into  engineering  units.  The  develt^iment 
of  the  stage  characteristics  used  in  this  rqiott  was  accomplished  by  Copenhaver  [2]  and 
O’Brien  and  Boyer  [14].  The  conversitm  of  high-response  data  into  engineering  units 
is  covered  in  Section  5. 

Measurement  Uncertainty 

The  method  used  to  determine  measurement  uncertainty  for  the  high-response  and 
time-averaged  acquisition  methods  was  defined  by  Abemethy  and  Thompson  [22].  The 
bias  and  precisitm  errors  were  cmnbined  using  the  equation: 

where: 

=  bias  oTors  based  on  calibration, 

ss  precision  errors  associated  with  the  transducer,  signal  conditioning, 
and  recording  devices, 

hs  =  stuttent  T  distribution  weighting  factor  based  on  the  number  of  sanqiles 
used  to  calculate  a  measurement  average. 

The  time-averaging  of  30  samples  resulted  in  t„  being  set  equal  to  2  for  the  time- 
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averaged  measoreinents.  For  this  report,  dw  weighting  factor  used  for  high-response 
measurements  was  also  2  because  each  measurement  was  an  average  of  500  samples. 
The  estimated  bias  errors,  precision  errors  and  measurement  uncertainties  for  the  time- 
averaged  meaairetnents  used  in  this  report  are  listed  in  Table  B.l.  The  esdrruued 
precision  errors  associated  with  the  raw  high-response  measurements  are  shown  m  Table 
B.2.  The  bias  errors  associated  with  the  high-response  measurements  could  not  be 
properly  estimated,  since  all  of  the  signal  conditioning  and  digitization  equipment  added 
an  unknown  offset  to  all  of  these  measurements.  However,  the  bias  errors  are  not 
important  here,  because  they  are  absorbed  during  calibratirm.  So  long  as  the  bias  errors 
are  assumed  to  be  constant,  the  tme  voltage  measurements  provided  by  the  high-response 
instrumentation  are  not  inqwttant. 


Table  B.1  Time-averaged  measurement  uncertainty 


Locatka  &  Type 

Range 

Bias 

Errer 

Predaieo 

Error 

MeasarenMot 

Uncertainty 

Inlet  static 
piessiiie 

(MW  iriMdicw  «1) 

34JkPa 

(3|ie4 

0.1329  kPa 

(0J)2217  pea) 

0.03851  kPa 

(pjOQ3S86pei) 

0.1712  kPa 
(002483  pei) 

Stage  10  inlet 
static  piessuie 

(new  inMdacw  #2) 

345  kPa 

(SOpei) 

1.275  kPa 

(0.1849  pei) 

0.2866  kPa 
(Or>4137  pei) 

0.1398  kPa 

(002027  pci) 

Discharge  total 
piessuie 

(new  nenedneen 
«3  AtS) 

690  kPa 

(iOOpei) 

2J48kPa 

(03695  pei) 

0.5715  kPa 

(008290  pei) 

2.793  kPa 

(04030  pei) 

Discharge  static 
pleasure 

(new  Beneitncwe 
#4  *#6) 

690  kPa 

(IOOpei) 

2.920  kPa 

(0/4336  pei) 

0.8956  kPa 
(0.1299  pei) 

3.426  kPa 
(04969  pei) 

I 
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Table  B.2  High-response  precision  errors 


Range 

(Wtor  •Ipeel 

Precision 

Error 

Probe  Transducer  #1 

5  volts 

7.26  mv 

Probe  Transducer  #2 

5  volts 

4.33  mv 

Probe  Transducer  #3 

5  volts 

7.S>4  mv 

Probe  Transducer  #4 

5  volts 

4.78  mv 

Probe  Transducer  #5 

5  volts 

13.25  mv 

Probe  Transducer  #6 

5  volts 

5.23  mv 
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APPENDIX  C:  Surge  Data 


The  high-response  surge  data  used  in  this  rqx>rt  are  presented  in  diis  appendix. 
Seven  plots  showing  pressure  and  discharge  valve  area  as  a  function  of  time  are  presented 
for  each  of  die  12  surge  events.  For  those  who  oppose  adoption  of  the  Systime 
International  d’Unit6s,  1  psi  »  6.895  kPa. 
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TabtoC.1  Surge  events 


PtonuRi 

Vokam 

Surge 

mMnOm 

SM 

NUmMT 

81%  * 

100% 

1 

17 

2 

18 

3 

19 

4 

20 

5 

70 

120% 

6 

28 

80% 

7 

90 

82% 

100% 

8 

7 

9 

15 

80% 

10 

87 

87% 

100% 

11 

8 

80% 

12 

86 

*  Stall  nomber  MicMM  (he  occataace  of  iIib  tar|e  or  rotatins  Maii  relative  to  oilier  tech  evenla  daring  the  entire  compieaaor 
tort.  For  exanpie,  rtall  namber  M  wae  (he  Mth  "rtaU"  of  te  toat. 


127 


SURGE  EVENT  «1 

(Stall  #17) 
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Figure  C.1  Surge  Event  #1  Data 
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SURGE  EVENT  #2 

(Stall  #18) 
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Figure  C.3  Surge  Event  #2  Data 
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Figure  C.4  Surge  Event  #2  Data 


SURGE  EVENT  «3 

(Stall  #19) 
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Figure  C.5  Surge  Event  #3  Data 


SURGE  EVENT  «4 
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Figure  C.7  Surge  Event  #4  Data 
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Figure  C.8  Surge  Event  #4  Data 
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Figure  C.10  Surge  Event  #5  Data 


SURGE  EVENT  #6 

(Stall  #28) 
Figures  C.11  andC.12 
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Figure  C.11  Surge  Event  #6  Data 


Figure  C.12  Surge  Event  #€  Data 


SURGE  EVENT  #7 

(Stall  #90) 
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Figure  C.13  Surge  Event  #7  Data 
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Figure  C.14  Surge  Event  #7  Data 


SURGE  EVENT  #8 

(Stall  #7) 


slj 

§■8  i 

«  I  8 

1^3 

Mi 

Ll.  00  o 


9  18  8  8  1 


142 


Figure  C.15  Surge  Event  #8  Data 
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ireC.16  Surge  Event  #8  Data 


SURGE  EVENT  #9 

(Stall  #15) 
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igureC.17  Surge  Event  #9  Data 
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Figure  C.18  Surge  Event  #9  Data 


SURGE  EVENT  «10 

(Stall  #87) 
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Figure  C.19  Surge  Event  #10  Data 
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Figure  C.20  Surge  Event  #10  Data 


SURGE  EVENT  #11 

(Stall  #8) 


Figure  C.21  Surge  Event  #1 1  Data 
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Figure  C.22  Surge  Event  #1 1  Data 
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Figure  C.23  Surge  Event  #12  Data 
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Figure  C.24  Surge  Event  #12  Data 


